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1, INTRODUCTION 


Soon after the discovery of the optical activity of quartz by Arago in 1811, 
Biot showed that the magnitude of the rotation in the visible region of the: 
spectrum is inversely proportional to the square of the wavelength, i.e., 
p=k/d*. The rotary dispersion of quartz has been, subsequently, the 
subject of numerous investigations. These have been reviewed by Sosman 
in his book on the properties of silica (1927) and by Lowry in his treatise on 
optical rotatory power (1935). More recently, the measurements have been 
extended into the ultra-violet upto 1525 A.U. by Servant (1938) and into 
the infra-red upto 9-74 by Gutowsky (1951). From these determinations 
it is evident that Biot’s Law fails completely as we move away from the 
visible region of the spectrum in either direction. Indeed, in the ultra-violet, 
the rotatory power increases from 6000 A.U. to 1500 A.U. nearly 40 times 
instead of 16 times as required by the inverse square formula. Several 
attempts have been made to represent the rotary dispersion numerically in 
terms of various assumed dispersion frequencies. A critical examination 
of these formule, however, shows that they are arbitrary in character and 
do not possess any real physical justificatior. 


A noteworthy fact about the optical rotatory power of quartz is that 
it increases appreciably. when the temperature of the crystal is raised, the 
rate of increase accelerating as the temperature approaches that at which 
the a-8 transformation in crystalline quartz takes place. At the transition 
point itself, there is an abrupt increment in the rctatory power. The tempera- 
ture-coefficient remains positive in 8-quartz but is very much smaller than 
in a-quartz. Though many experimental studies dealing with these effects 
are on record, no attempt appears to have been made to give a theoretical 
explanation of the influence of temperature on the optical activity of quartz. 


It is proposed to discuss in this paper, firstly, the manner in which the 
optical rotatory power of quartz varies with wave-length and secondly, its 
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variation with the temperature of the crystal. It emerges from the discussion 
that the two topics are closely interrelated. 


Other crystals, viz., cinnabar and sodium chlorate, are known which 
are optically active and in which also the temperature coefficient of the 
optical rotatory power is positive. These cases will be discussed in later 
papers. 

2. THE EARLIER FORMULA FOR THE ROTARY DISPERSION 


The first attempt at a theoretical formula to fit the rotary dispersion 
of quartz was made by Drude in his ‘* Optics ” where he proposed the form 


P “es ee (1) 


where p is the optical rotatory power (in degrees per millimetre), Q, is a 
constant corresponding to the characteristic absorption wavelength 4, 
(in microns), the summation being performed over all the r absorption 
frequencies. Drude found that the then available data could be fitted by 
a two-term formula, where 


Q; = 12-2 Q, = — 5:046 
A,?= 0-010627 A,? = 
the negative term being unambiguously indicated to be necessary. 


Since then several other modified formule have been suggested by 
various authors. Lowry (1912) noticed that Drude’s formula is not 
accurate enough to fit his first series of extended measurements and hence 
altered it to 

11-6064 4-3685 . 13-42 
P= = 0-010627 FO 78-22 2) 
the additional third term corresponding to the infra-red absorption band. 
His later still more accurate readings upto 2263 A.U. however required still 
further modification and his final formula is 


9-5639 23113 
EE Tn 
P = y3—0-0127493 — X®— 0-000974 — 9° 1905 (3) 
where instead of the single absorption band at 1031 A.U. occurring in the 
previous formule, two new characteristic wave-lengths 1129 A.U. and 
312 A.U. have been used, the effect of the infra-red being substituted by 


a constant. 


Neither of these latter formule, however, represent the results of 
Duclaux and Jeantet (1926) which extend upto 1854A.U. and therefore 
another formula was worked out by Bradshaw and Livens (1929), viz., 
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_ 845-694 040235, 
? = (X2— 0-01274912) ~~ (A2@— 0-01274912)? 
_ 838-4320 (01331233, «4305794 
(A2— 0-0120800) ~ (A2— 0-0120800)2 ~ (A2— 80-00) 
2119-117 | 
> (= 80-0002" 7 


It will be noticed that in this formula, the optical rotatory power is expressed 
as a difference of two very large terms of nearly equal magnitude but of 
opposite signs involving two close absorption wave-lengths respectively at 
1129 A.U. and 1099 A.U., and that the influence of the infra-red term has 
been retained. Though some additional terms of a different type occur 
in this formula, it is clear that the terms of the Drude type appearing in it 
are the priacipal ones. The formula. however, gives values which deviate 
notably from the newer extended measurements of Servant (/oc. cit.) in 
the remote ultra-violet. Radhakrishnan (1947) has therefore fitted these 
results with still another rotary dispersion formula, viz., 


_ 4948 4617 
° = = 0-014161 ~ = 0-011236 ~ A*=0-000974 
0-000815 
~ (R= O-o195;" — 9° 190° ©) 


This formula is largely a refinement on the work of Lowry and Coode Adams. 
The absorption at 312 A.U. is still retained but the single term at 1129 A.U. 
has been split into two terms at 1060 A.U. and 1190A.U. In addition, a 
small negative term at about 1400 A.U., almost at the very limit of conti- 
nuous absorption in quartz has been introduced. This formula fits the entire 
data fairly well. 


Thus, we see that a variety of formule involving numerous hypothetical 
absorption frequencies have been suggested to numerically represent the 
rotary dispersion of quartz. On examining these formule, we notice that, 
without exception, all the authors have been obliged to introduce negative 
terms involving absorption wave-lengths lying in the ultra-violet, in order 
to fit the data over the whole spectrum. No physical justification has been 
given for this procedure. The necessity for assigning a negative strength 
to a term becomes clear when we consider a single-term formula of the type 
(1) as suggested by Drude. When A, = 0, this reduces to Biot’s Law. When 
A, takes finite values, the rotatory power varies faster with wave-length than 
required by Biot’s Law. The rotatory power has been calculated for several 
wave-lengths giving various values to A, in this formula and shown in 
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Table I. These values of A, have been chosen at random over the entire 
region of continuous absorption in quartz which extends below 1436 A.U. 
(Tousey, 1940). In each column, the rotatory power has been fitted with 
the experimental value for A= 0-6708u. (Throughout, p is expressed in 
degrees per millimetre.) 

TABLE | 


Rotary dispersion using one-term formula of Drude type 


p (Calculated) 








r (E ae | 
(in p) mee ee Siti 4 a. ae Oe ae — 
| , | mental) | 
hy =0 Ay = 03 | Ay = -06u | Ay ="10u Ay= "12 
| a. 4 
1525 220 332 75-6 549 813 776 
1625 281-8 29} -2 23+7 443-5 600-0 : 599°5 
1825 293-4 229-1 248-5 312-2 381-0 | 391-5 
2359 133-7 135-6 141-8 159-4 174-6 =| 180-4 
-3728 58-6 | 58-8 54-5 56-4 57-8 i 
4678 34-0 | 4-1 34-3 34°8 35-2 35-6 
- 6708 16-54 | 16-54 16-54 16-54 16-54 | 16-54 

















[It is seen 1 that for ail values of r, ‘except 1, = 1200 A.U., ‘the calculated 
rotatory power is much lower than the observed value over the whole range 
of wave-lengths. The same is true even for A,= 1200 A.U.' except in the 
remote ultra-violet region where the calculations exceed the observed data. 
On a careful study of this table it can easily be understood why no combina- 
tion of positive terms can ever accurately represent the rotary dispersion 
of quartz. For, even if at the two extremes of the spectrum the values be 
accurately fitted by a suitable choice of terms, there will be a large region 
in the visible and near ultra-viclet for which the calculated values will be 
less than the experimental values. This difficulty has been overcome by the 
previous authors by expressing the rotatory power as a difference of two 
terms. Evidently, these formule are purely empirical in character and do 
not give a real explanation of the facts. Hence, we may conclude that a 


rotary dispersion formula of the type suggested by Drude is not appropriate 
in the case of quartz. 


3. A New FORMULA 


The present author has found that the rotary dispersion of quartz is 
accurately represented throughout the visible and ultra-violet range of the 


spectrum by a simple formula containing only two constants, viz., 
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ki? 
P= (1,5 (6) 
It is seen that the formula approximates to Biot’s Law when A >> Aj. 
Also for a given value of A,, the rate of variation of the rotation with wave- 
length from this formula is much greater than from the Drude formula. 
The constants which fit the experimental data are 


k = 7-186 and A,= 00926283 p. 


The rotatory power calculated from this formula are given in comparison 
with the experimental values in Table II. 


TABLE II 
Rotary dispersion using author’s formula 














A (in w) p (experimental )| p (calculated) X (in pw) |p(experimental Ne (calculated) 
J | | | 

+1625 776 7176 + 291216 104°97 104-88 
+1575 7124 724 - 340365 72-46 72-35 
+1625 599-5 597-2 *372762 58-84 58-75 
*1675 533-0 531+5 + 435834 41-55 41-50 
+1700 504°5 503-0 -467816 35-61 35°57 
+1750 453-5 452-8 -479991 33°68 33-65 
*1775 431-5 430-7 - 508582 29-73 29-72 
-1800 410-5 410-3 -520908 | 28-25 28+24 
+1825 391-5 391-5 +§46074 } 25-54 25-55 
+1854 370-9 371-3 -579066 22-55 22-57 
+1935 322-76 322-89 *643847 18-02 18-08 
- 198979 295-65 295-82 - 670785 16-54 16-60 
+217402 226-91 226-96 -810 | 1l-1l 11-24 
+ 226909 200-90 200-96 -940 8-14 8-29 
+ 235923 180-43 180-44 1-170 5°18 5-32 
+ 241331 169-68 169-70 1-342 3°89 4-03 
+ 262830 135-66 | 185-62 2-500 0-98 1-15 
*273955 122-12 122-06 3-210 | 0-52 0-70 





While the previous formule, with the single exception of Radha- 
krishnan’s fail in the remote ultra-violet, this one fits very well from 6000 A.U. 
to 1525A.U. In fact from 5800 A.U. to 1800 A.U., the maximum error 
is only about one part in a thousand, being sometimes positive, sometimes 
negative. Between 1800 A.U. and 1525A.U., the deviation is slightly 
larger for certain wavelengths, but even so it is not of a progressive nature. 
There is a small discrepancy as we proceed from the red end of the spectrum 
to the infra-red, the error uniformly increasing with increase of wave-length. 
This is no doubt due to the contribution of the infra-red absorption bands, 
which have not been included in the above formula. As a matter of fact, 
Gutowsky (Joc. cit.) has found that quartz exhibits anomalous rotary dis- 
persion beyond 3-7 u, the rotation, instead of decreasing as is to be expected, 
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increasing nearly four times from 3-7 to 9-74, thereby indicating that 
the presence of absorption bands in the infra-red has an influence on the 
optical rotatory power. 


It can easily be shown that the single characteristic frequency used in 
the above rotary dispersion formula lies within the region of the ultra-violet 
spectrum in which the absorption frequencies of quartz appear. A one- 
term dispersion formula for the ordinary refractive index has been fitted 
up using this same frequency. The formula proposed is of the Sellmeyer- 
Drude type which is known to vield satisfactory results in the case of quartz 
(Radhakrishnan, 1951). The formula is 





‘ 1-35 A? 
Ny? — 1 = vs-—aao 6983)" (7) 
A?— (0-0926283) 
The calculated values are given in Table III. 

TABLE III 
Dispersion of quartz (for explanation see text) 
» (in #) 2 (experimental) (calculated) | 
-1854 1-676 1-673 | 
| 
“2413 1-605 1-607 | 
-3404 1-5675 1-5678 
5086 1-548 1-548 | 
- 6708 1-541 1-541 | 
- 8000 1-538 1-539 | 


1-2000 1 +532 1-536 | 
| 


We see that the error in the calculated value over the whole range is only in 
the third decimal place. This indicates that the characteristic wave-length 
of 926 A.U. is roughly the same as the effective average wave-length of the 
absorption spectrum of quartz in the ultra-violet. 


4. THE TEMPERATURE VARIATION OF THE OPTICAL ROTATORY POWER 


The rotary dispersion formula (6) may be expressed in the form 


9 
av~ 


~ Oe 7 


p 


+ 


where ‘a’ is a constant such that k = ad,‘ and », is the characteristic 
absorption expressed in wave-numbers. The optical rotatory power of 
quartz is known to increase with rise of temperature even when the expansion 
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along the optic axis has been corrected for. To calculate this increment 
let us differentiate equation (8) with respect to temperature. We shall 
simply consider ‘a’ to be constant with temperature, but », to vary with 


0 
it, so that 
= dav*y 
~ (% ets Xo (9) 
where 
» 1 dry __ 1 da, 
ee aa ra 
Hence the temperature coefficient of the optical rotatory power 
ldp _ — 42? 
p at ~ (2 —1,”) Xo- (10) 


Molby (1910) has made precision measurements of the optical rotatory 
power of quartz at temperatures of — 118°C. and 20°C. from 4000 A.U. 
to 6700A.U. With the aid of equation (9) and the rotatory dispersion 
formula (6), the values of do/dt have been theoretically accounted for over the 
whole range of wave-lengths. The value of x used is 2-7 x 10-5. The 
theoretical values are given in Table IV. 


TABLE IV 


Temperature variation of optical rotatory power 











| p (experimental ) | dp 
A dp ~ (theor.) 
| (in p) di (expt. ) dt a 
20°C, | —188° C. | x=2e7x10™ 
| 
| ; 
| 
+4047 48-842 | 47-675 5-6x 107% 56x 1073 
+4359 41-441 40-469 | 4°7xX10-% 4°7x1078 
+4670 35-529 34-696 4-0x 1073 4-0x 1073 
+4800 33-593 32-818 3-7x 107-3 3°8 x 1073 
| -5086 29-656 28-983 3-2x 1073 3-3 x 1073 
-6461 | 25-472 24-885 2-8x10-3 2-8 x 1073 
| +5893 21-686 21-188 2-4x10-° 2-4x 1073 
* 6438 18-003 17-603 1-9x1073 2-0x 1073 
*6708 16-493 16-145 1-7x1073 1-8x1lo 3 





| 


ae Bn | 


Eicon pconeniacicadad elec acacia atic ate berlin 


The agreement between the theoretical and experimental values 3 is satisfactory. 
From equation (10) we notice that the theory predicts that the temperature 
coefficient of the rotatory power should exhibit an increase with decrease 
of wave-lengths. This in fact is known to be true from Molby’s measure- 
ments and also from the observations of Soret and Sarasin (1882). 





Le Chatelier (1889) has made measurements of the rotatory power over 
wide ranges of temperatures both below and above the «-f transformation 
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point in quartz. Using a value of y = 4-5 «% 10-5, dp/dt has been theoretically 
calculated as before. These are given in Table V. 
TABLE V 


Tenmiperature variation of optical rotatory power 


| 
| 


A (in 4) “ (experimental) . (theor.) 
x=4-5x10-° 
+279 27-5 x 107* (20° — 475° C.) 24 x 107" | 
+448 7-0 x 107° (z0°—415° C.) | 7*3x 107° | 
-500 5*4x 107" ( m ) 5°7X10-" | 
“518 5-0 x 107" ( re ) 5°3 x 107° 
| 589 4-3 x 107" ( : ) | 4-0x 1073 
- 656 3-4 x 10-*( m ) |  8-2x10-2 





The value of x in this case has increased as the temperature is higher. In 
fact, as will be seen later, x is found to have almost the same value and 
increase in roughly the same proportion with temperature, from the refrac- 
tion data as well. Thus the increase of temperature coefficient with increase 
of temperature and with decrease of wave-length have both become intelligible 
from these calculations. 


5. RELATION TO ORDINARY DISPERSION AND ITS THERMAL VARIATION 


It would be of interest te see how the results obtained in the previous 
section are related to ordinary dispersion and its variation with temperature. 
Dispersion formule for quartz have been suggested by manv authors in the 
past. We shall consider nere only that put forward by Radhakrishnan 
(loc. cit.) which has the merit of fitting both the ordinary and extraordinary 
indices very accurately. We shall quote only that for the ordinary index, viz., 

0- 663044 A? 0-517852 A* 0+ 175912 r2 
2 10-0600)? © A2— (0-1060)2 © A2— (0-1190)2 
0-565380.42 | 1-675299 A? an 
A?— (8-844)* * A?— (20-742)*° ; 


M,*— 1 = 


In this formula, the absorption at 600 A.U. is a hypcthetical one, while 
those at 1060 A.U. and 1190 A.U. correspond with the reflection maxima 
observed by Tousey (/oc. cit.). In addition, there are two small infra-red 
terms which serve to give a very accurate fit over the entire range of available 
data, 
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It is quite possible, of course, that the characteristic wave-length at 
9296 A.U. which determines the optical rotatory power also plays a sensible 
role in the ordinary dispersion of quartz. Accordingly, a dispersion formula 
for the ordinary index of quartz has been fitted up in which this wave-length 
has been included. The absorptions at 600, 1060 and 1190 A.U. have been 
retained, but an additional one at 926A.U. has been introduced. The 
formula 1s 


2) — 0765042 010000384 
Mo * = 2 (-0600)2 * AZ— (0-0926283)2 ~ A2 — (0- 1060)? 
+2228 A2 


. , — ° 2 e 

12 — 10-1190)2 0-01193 A?, (12) 
The agreement between the calculated and the experimental values can be 
seen in Table VI below. 


TABLE VI 


Dispersion of quartz ( for explanation see text) 








d (in w) ny (expt. ) | my (calc.) | 
- 185467 1-67578 1+ 67578 
+ 214439 1-63039 1-63038 
+ 2912358 1-58098 1-58096 
- 340365 1-56747 1-56747 
-479991 1 -55012 1-55012 
+ 58929 1+54425 1+54425 
- 670786 1-54147 1-54146 


The accuracy of this formula is seen to be as good as any of the formule 
suggested previously. 


Measurements of dn/dt have been made by Radhakrishnan (/oc. cit.) over 
a wide range of temperatures for several wave-lengths. Using the data for 
the ordinary index, the y’s for the several dispersion frequencies occurring 
in Formule (11) and (12) have been calculated in the manner suggested by 
Ramachandran (1947). Since the contribution of the infra-red frequencies 
is small, their variations with temperature have been neglected. It is found 
that the data agree well when it is assumed that the extreme ultra-violet 
term at 600 A.U. is unaffected by temperature in both the formule. In 
Formula (11) we shall assume that the y’s for 1060 A.U. and 1190 ALU. 
are the same, as these two terms are very close to each other, Similarly, 
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in Formula (12) we shall assume that the x’s for 926, 1060 and 1190 are the 
same. The ’s calculated from the two formule are tabulated below. 


TABLE VII 


Temperature variation of refractive index 











) 
| 
Temperature 50° C. Temperature 400° C. 
‘ | 7 = =e “eames 
(in w) _ 2 195 xx105 | yxloe | a nat x x105 xxl | 
at Formala (11)|Formula (12)! at Formula (11)|Formula (12)! 
| (expt.) | (expt. ) 
| | aw es oe 
+2537 0-29 | 2-3 2-3 1-00 3-5 3-4 
| 
+3650 0-54 | 2-2 21 1-32 393 3+3 
-5461 0-62 2+] 2-1 1-44 3°3 3-2 
| +5893 | 0-63 2-1 2-1 1-46 | 3-2 3-2 


| = oe = 








It will be noticed that the values of x obtained from formule (11) and (12) 
are almost the same. It will also be seen on comparing the above table 
with Tables TV and V that the y’s estimated from the data on the thermal 
variation of refraction and the thermal variation of rotation are about the 
same and also increase in the same manner with rise of temperature. 


6. DISCUSSION 


The fact that the rotary dispersion of quartz is expressible by a formula 

of the type 
P (A?—A, 22 

indicates that its optical activity is a result of an interaction or coupling 
between the successive polarisable units in the structure of the crystal 
excited to similar upper energy levels, or in other words, an interaction in 
the nature of resonance. In their article in the Handbuch der Physik, Born 
and Goppert-Mayer (1933) arrive at the general formula 


aad Pw , @ . » FY . G 


; (A2—A?)? A2—d,? ss i : (A2?— 2,7)? \2— r,2” 
It will be noticed in this formula, that in addition to terms of the Drude type, 
certain quadratic terms, of the same form as suggested in this paper, are 
also present. In view of the close agreement between the experimental 
values and the values calculated from the author’s formula (see Table II), 
it appears improbable that terms of the Drude type play any observable role 








) 


1 
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in determining the optical rotatory power of quartz in the region of the visible 
and the ultra-violet part of the spectrum. 


In conclusion, | express my sincere gratitude to Prof. Sir C. V. Raman, 
F.R.S., N.L., for the kind advice and the encouragement he gave me during 
the course of this investigation. 


7. SUMMARY 


N@nerous formule have been proposed by several authors to account 
for the rapid increase in the rotatory power of quartz with decrease of wave- 
length. These formule generally consist of a large number of terms both 
positive and negative, principally of the Drude type, involving various 
hypothetical characteristic frequencies. An examination of these formule 
shows that they do not possess physical justification and are arbitrary in 
character. In the present investigation it is shown that the data are well 
represented by the single term formula 

kA? av? 
(2—A,2)2 or (v,2— v2)” 
where k = 7-186 and A, = 0-0926283 4. Assuming that ‘a’ is a constant 
with temperature and »v, to vary with it, the thermal variation of the optical 
rotatory power has been quantitatively accounted for, vy, being found to 
vary in the same manner with temperature as estimated from the thermal 
variation of the refractive indices of quartz. 


o> 
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in Presence of Finely Divided Copper 


By G. V. JADHAV AND V. N. DESHPANDE as 


(Karnatak College, Dharwar, and Organic Chemistry Laboratory, The Institute of Science, Bombay) 


Received November 10, 1951 
(Communicated by Dr. Mata Prasad, F.A.Sc.) 


Ir was observed by Hirwe, Jadhav and Chakradev! that the reaction of 
thionyl chloride and methyl §-resorcylate in presence of finely divided copper 
was very vigorous and no definite product could be isolated. It has, how- 
ever, been now found that diary! disulphide (dithio-ether) could be isolated 
when the reaction was carried out at a lower temperature. 

The point of linking of the two nuclei in this disulphide was arrived 
at by treating it with nitric acid in presence of sulphuric acid below —5S° 
when 5-nitro-8-resorcylic acid was obtained, which showed that the two 
nuclei were linked through sulphur atom which occupied the position para 
to —OH group and meta to —COOCH, group. If, however, this nitric acid 


reaction was carried out at room temperature 3: 5-dinitro-B-resorcylic acid 
was obtained. 


The disulphide should have the constitution as ns = § (where R 
is the nucleus) and not R—S—S—R as it gave 2: 2’-dimethoxy-4: 4’-di- 
hydroxy-5: 5’-dicarbmethoxy-diphenyl thioether when methylated with 
dimethyl sulphate in presence of alkali. This thioether gave 2: 2’-dimethoxy- 
4: 4’-dihydroxy-5: 5’-dicarboxy-diphenyl thioether identical with the one 
got by Hirwe, Jadhav and Chakradev (Joc. cit.). If the thioether would 
be supposed to have the constitution R-—-S—S—R both sulphur atoms 
should have been lost simultaneously. 

As the same disulphide (dithioether) was formed by the acticn of sulphur 
monochloride as wel! as that of sulphur dichloride, the mechanism of the 
reaction may be explained by considering that thionyl chloride first decom- 
poses as (1)? or (2) 

3 Cu + 4 SOCI, —3 CuCl, + 2 SO, + S,Cl, (1) 
Cu + 2 SOC], CuCl, + SO,-+ SCI, (2) 
and then 2 SCI, — S.Cl. + Clo. 
114 
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As the yield of the disulphide with sulphur dichloride and sulphur mono- 
chloride is poor as compared with that obtained with thionyl chloride and 
copper, copper may be acting as catalyst also in this reaction. Further as 
the yield with sulphur dichloride is better than with monochloride, chlorine 
formed by its decomposition may be acting as a catalyst. 

Sulphur monochloride has been known to form dithioethers (disulphide) 
by Naik*® and Naik and Bhat.‘ 

The diacetyl and dibenzoyl derivatives of the disulphide were also 
prepared. 

EXPERIMENTAL 

2:2'-4: 4'-Tetrahydroxy-5: 5'-dicarbmethoxy-diphenyl disulphide (I).— 
Methyl-f-resorcylate (6 g.) and thionyl] chloride (25 c.c.) were mixed together 
and the mixture was surrounded by ice. Finely divided copper (8 g.) was 
gradually added to it. Fumes of hydrogen chloride and sulphur dioxide 
were evolved. The reaction mixture was protected from moisture and kept 
overnight in refrigerator. It was then mixed with dry chloroform and heated 
on a boiling water-bath for two hours and the hot solution filtered. A pasty 
mass was obtained after the removal of the liquid, which turned into solid 
after repeated washings with acetone. The solid was then boiled with 
carbon disulphide to remove free sulphur and finally crystallised from a 
mixture of nitrobenzene and acetic acid, m.p. 245-6°. Yield: 4g. Found: 
S, 16-1; C,,H,4O,S2 requires S, 16-08 per cent. 

5-Nitro-B-resorcylic acid.—The disulphide (I) (0:5 g.) was dissolved in 
sulphuric acid (5c.c.) and the solution was surrounded by freezing mixture. 
Nitric acid (con.) (5¢.c.) was gradually added to it, not allowing the tempe- 
rature to rise above —5°. The reaction mixture was kept in ice after the 
addition was over and then diluted with water when solid was obtained, 
which crystallised from hot water. m.p. 215-6°. This showed no lowering 
in melting point when mixed with a genuine specimen. 

3: 5-Dinitro-B-resorcylic acid.—The disulphide (I) was dissolved in 
sulphuric acid (5c.c.) and nitric acid (con.) (5c.c.) was added to it when a 
vigorous reaction commenced which was allowed to proceed for ten minutes. 
The precipitate obtained on diluting the mixture with water was crystallised 
from dilute acetic acid. It sintered at 190° and melted completely at 204° 


(decomposition). It showed no lowering in melting point when mixed with 
a genuine specimen. 


2: 2'-Diacetoxy-4: 4'-dihydroxy-5: 5'-dicarbmethoxy-diphenyl disulphide 
(JI).—It was prepared by acetylating (1) in presence of pyridine. Crystallised 
from glacial acetic acid, m.p. 213-4° (decomposition). Found: S, 13-4; 
CoH; s019S2 requires S, 13-3 per cent. 
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2: 2'-Dibenzoyloxy-4 : 4'-dihydroxy-5 : 5’-dicarbmethoxy-diphenyl disulphide 
(11T).—It was crystallised from acetic acid. It shrinks at 170° and melts at 
191-2° (decomposition). Found: S, 10-3: Cy gHo:O0, 9S. requires S, 10-6 
per cent. 

2: 2’-Dimethoxy-4 : 4'-dihydroxy-5 : 5'-dicarbmethoxy-diphenyl — sulphide 
(1V).—This was obtained when the disulphide (1) (0-5 g.) was heated on a 
boiling water-bath with dimethyl sulphate (5 c.c.) in presence of 20% solution 
of potassium hydroxide (40c.c.). It crystallised from acetic acid in colour- 
less needles which did not melt upto 300°. Found: S, 7:7; C,,H;,0,S 
requires S, 8-1 per cent. 

This sulphide (IV) was hydrolysed to 2: 2’-dimethoxy-4: 4’-dihydroxy- 
5: 5’-dicarboxy-diphenyl sulphide by boiling it for 2 hours with 10 per cent. 
sodium hydroxide solutien. It crystallised from acetic acid, m.p. 224—-5°. 
Found: S, 8-6; C,,H,,O,S requires S, 8-7 per cent. Hirwe, Jadhav and 
Chakradev (/oc. cit.) gave m.p. 226°. 

It showed no lowering in melting point when mixed with the genuine 
specimen. 

The same disulphide (TI) was also prepared by allowing sulphur dichloride 
(5c.c.) to interact with methyl B-resorcylate (4 g.) in presence of dry ether 
(20 c.c.) and allowing the reaction to proceed for two days at room tempe- 
rature. The semi-solid obtained after the removal of ether was stirred with 
carbon disulphide when grey solid separated. It crystallised from acetic 
acid, m.p. 245-6°. Yield, 1:7g. Found: S, 16°03; C,,H,,0,S. requires 
S, 16-08 per cent. 

Sulphur monochloride also gave the same product when the reaction 
was allowed to proceed at room temperature for two days, but the yield was 
very poor (0-7 g. only from 4g. of ester). 


SUMMARY 


The disulphide was obtained in this reaction and it was found to be 
unstable towards alkali. Its characteristic derivatives such as di-acetoxy 
and di-benzoyloxy were prepared. 
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1. INTRODUCTION 


TITANIUM DIOXIDF exists in three polymorphic varieties, namely, rutile and 
anatase (tetragonal) and brookite (orthorhombic). These structures are 
alternative forms of 6:3 co-ordination. All the forms are strongly absorb- 
ing below the blue region, the absorption being usually enhanced by 
impurities. The refractive indices are of the order of 2-5 to 3-0, with 
correspondingly large values of dispersion and birefringence. Measure- 
ments of dispersion have been made by many workers chiefly Schréder 
(1928) whose data have been used in the International Critical Tables 
(Vol. VIN) to furnish dispersion data and estimates of dn/dt. The principal 
thermal expansion coefficients of rutile and anatase have been measured 
by Fizeau (1866) and those of brookite by Schrauff (1884). These data 
on refraction and expansion have been used here to compute dispersion 
formule for the three forms and calculate the temperature shift of the 
ultra-violet dispersion frequencies. Finally, the refractive indices of rutile 
have been calculated by Bragg’s method and good accord obtained with 
the observed values. 


2. DISPERSION FORMUL2 FOR TITANIUM DIOXIDE 


Schréder’s dispersion data (Schréder, Joc. cit.) range from 4,500 to 
7,000 A and are sufficiently extensive to locate the nearest ultra-violet dis- 
persion frequency and suggest a correction for infra-red effects. Thus the 
proposed formule are of the form 

. 
n? —— 1 =24-r A2— dj? un cA?. (1) 
The value of Ay can be made the same in all the three modifications, namely, 
0-2870. The other constants of the dispersion formule and the agreement 
between calculated and observed indices are shown in the succeeding 
tables. 
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TABLE | 


Dispersion Formule for Anatase 























Index | a / 
w 2 +4485 2+3998 | 0-10 
‘ 2+8225 1+8492 0+18 . 
| \ (in w) n (calc. ) 1 (obs. ) | z (calc.) 2. (obs.) 
| 0+ 43583 2+ 7688 2.7688 2.6576 2+6576 
0-49161 2-6584 2+6579 2+5688 2+ 5679 
| 0-54607 2-5948 2-5951 2-5161 2-5166 
0-57907 2.5679 2.5683 2-4936 2-493 
0+60727 2+5498 +5498 2-4793 24793 
0-67163 265179 2-5179 2°4518 264507 
0+ 70820 2-5072 2-5039 2+4404 
TABLE Il 
Dispersion Formule for Rutile 
| 
Index a | b 
ARETE RATE SET ADE, LO ORT 
| 
w 2+ 2868 207275 0-10 | 
P 2+6034 3+7199 0-20 | 
| 
| x (in «) | a, (cale.) | n _, (obs.) | n. (calc.) n, (obs.) 
a | 
0)+54607 2-6505 ©. 6505 2-9467 | 29467 
0+57907 2.6211 2-6214 2- 9085 | 2.9100 
| 0-60727 2-6001 2+ 6003 2-9842 | 2.8842 
| 0-6 7163 2-5643 2.5648 2.8307 | 2.8397 
0+69075 22-5555 25555 2.8294 2.8994 
0- 70820 25495 25495 








net 


Si reer 








eS NTI 
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TABLE III 


Dispersion Formula for Brookite 





Index a b é 

a ..| 269858 2+1036 0-18 

| 8 ..| 268057 2-8130 0-16 

rr ve) 24285 2-9641 | 0-05 
\ (in pw) ne (calc. ) na (obs.) mg (calc. ) Mt obs.) | % (calc.) | tt, (obs.) | 
\ ‘aeeaiak: teciaiaaee anand i 5 eumaaecs 
043583 2+ 7686 2+ 7695 2+ 7838 2+7837 2-9419 2-9416 | 
0-49161 26708 * 2-6717 2-6770 2-6769 2-8046 | 
0+54607 2.6150 2-6154 2.6159 2-6160 2+7402 2+7402 | 
057907 2.5910 2-5904 2+5895 2+5900 2.7094 | 2.7091 | 
060727 25743 2+5739 2-5718 2 +5718 2-6882 62-6882 
0+67163 25441 |  2+8443 2+5405 2+5403 2-6518 | 2-6519 | 
0+ 69075 2+5369 2+5376 25331 2-5328 2 +6433 2-6429 | 
0-70727 2+5312 2+5317 2:5279 | 2+5265 | 





A perusal of these tables shows that the formule give results usually 
correct to the third decimal place, i.e., within the limits of experimental error: 
The nearness and the large strength of the term at 2870 A are responsible 
for the high refractive indices and dispersion. Since the constants a and b 
are approximately equal, the farther ultra-violet absorptions have together 
about the same strength in refraction as the nearest term. The infra-red 
correction term c has values of the order of 0-1 to 0:2, as compared with 
0-1 in quartz. Since the Raman spectrum (Narayanan, 1950) of rutile and 
the absorption data of titanium dioxide do not indicate any near infra-red 
absorption, the infra-red vibrations must have very large strengths. This 
is supported by the exceptionally high dielectric constants (of the order of 
50 to 100) of all the three forms at room temperature (Eucken and Buchner, 
1934). 


The near ultra-violet dispersion frequency is the same in rutile, anatase 
and brookite and hence the actual absorption peaks in this region must be 
nearly the same. This fact is understandable when we remember that in 
A2 
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all the three modifications, each oxygen has for its nearest neighbour three 
titanium atoms at the corners of a distorted equilateral triangle and each 
titanium atom has six oxygens as nearest neighbours at the corners of a 
distorted regular octahedron. The dispositions and distances of these 
atoms differ but little in the three forms (Bragg, Atomic Structure of 
Minerals, p. 104). Consequently the dispersion frequency is also nearly 
the same. It may be remarked that the situation in crystal and vitreous 
quartz is analogous. For both these, Tousey (1940) —— peaks of 
selective reflection in the same wavelength regions. 


3. THERMO-OpTiIC BEHAVIOUR OF TITANIUM DIOXIDE 


The thermo-optic coefficients dn/dt for the principal refractive indices 
of the three forms computed from Schréder’s data are given in‘the Inter- 
national Critical Tables. The values are extremely rough and can be used 
only to make approximate calculations on the basis of the theory of thermo- 
optic behaviour for anisotropic crystais (Ramachandran, 1947). 


Anatase.—In terms of Ramachandran’s theory, we can form for a uni- 
axial crystal the functions 


dn ‘ AP Dis 
Fu, = 2g GE + ¥% (Ma? 1) = Aw =) +E 7a? (2) 


dn, ‘ AP 2b, v9! 
= ote “ a = —_.—_-——_—. ie 

. n, at + 4 (n, 1) A.+ yyi— v2 | (¥% ae v2 
where y; is the coefficient of volume expansion of the crystal at the tempe- 
rature t, A,, A, are contributions to dn/dt from remote ultra-violet terms 
and AP represents a factor controlling the transfer of transition probability 
of the term at 2870 A from axial to perpendicular directions as the temperature 
is increased. The last term gives the contribution to dn/dt from the tempe- 
rature shift of the term with v, at 2870 A, the shift being defined by 

1 dy 1 da& 

as ie (3) 


Xy 


From equation (2) one can climinate AP and yx, in turn, obtaining the 
results 

2 (b+ 2By) Yo! 

iS + y 4 = A. + 2Aw + (v4? ail “aye : 


Pi tA A— Auk sae be. aP. 2(b) 


X, 2 (a) 


On plotting the left-hand side of Eq. 2 (a) against (v_2— v*)-®, one should 
obtain a straight line, the slope of which gives X, and the y- intercept 
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(A,+ 2A,,)- Similarly from Eq. 2(b) one can graphically estimate AP and 
(Ab, — Awb.). A, and A, can also be separately evaluated. The points 
show considerable scatter in these graphs, but from the mean straight lines, 
one obtains the values 


AP = 120 x 10°, A,= — 5:0 x 10°, A,= -- 15-0 x 10%, 
yy = 1:23 x 10-5. 
Using these constants, the calculated values of F,, and F, are compared, 
in Table IV with those directly obtained from dn/dt and y. The value of y 
was taken as 2°18 x 10-5 as obtained from Fizeau’s data extrapolated to 
90°C., the mean temperature of Schréder’s observations. 


TABLE IV 
Thermo-Optic Coefficients of Anatase at 90° C. 





| 


| A(in uw) | 10° Ge 108 ae 10°-F,, (obs.)|10°-F., (cale.) 10°-F, (obs.)| 10°+F (calc.)| 
0-4358 1-8 rO-4 4-6 4°07 15-4 15+4 | 
| 4916 264 —0-6 0-5 7 eS 8-9 | 
| 0+5461 se | -00 | ee ee 6-0 | 
| +6073 ~2+4 “153 os | -3 | 48 4 | 
| 0+6908 2.6 1-8 ~2+5 25 8| Bel 2-4 

| 


Table LV shows that the formule fit well the general trend of the 
observed values, but the casual errors in the latter are so large that a more 
accurate picture is not possible. It is for this reason that dn/dt has not been 
explicitly evaluated, the functions F,, and F, being more suitable for com- 
parison between theory and experiment since they are more smoothed out. 


Brookite.—\n this crystal, the refractive indices are very nearly the 
same and so too their temperature coefficients. One can therefore form 
the functions F,, F;, F, and write 


(F, + F;) ; AP _ (ba+ bs) vo*x,\ 
5 —* = 4(A, + As) + T(v,2— vi) 7 (v2 vty? | 
AP sod xy | (4) 


Fe A, aT Ge 


Using the value ygg= 5-57 x 10-° calculated from the observations of 
Schrauff, one obtains the coefficients 


S(A,+ A;) 0; A.= 10 x 10°; x,= 2-47 x 10°; AP = 440 x 10% 
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The calculated and observed values of (F,,+ F;)/2 and F, are shown in Table V. 


TABLE V 
Thermo-Optic Coefficients of Brookite at 90° C. 











| 5. 5. = | 
A(inw) | 10° | 108. 105. e (rat FN Fat Fore (Cal — | 
m | ee | 
0-4358 5+ | 5-3 66-3 66-3 aes | 
0-4916 3+7 | 3-0 —8-4 52-9 523 —11-4 | — 8-8 | 
| 0-5461 24 | 2 —9+1 45-9 45-3 | -12-2 | -13+6 | 
| o-s70 | 24 | 19 | 8-9 43-0 42-0 | -12-4 | -129 | 
0-678 =| 00 | ek | 89 38-0 39-1 145 | -1366 | 

| 0-7082 17 | kes ” 37-3 37-7 

| | 














As in the case of anatase, the calculated and experimental values have 
the same general trend. Owing to the large casual errors in the latter, a 
more accurate treatment is not possible. The values of x, and the transfer 
coefficient are twice the corresponding terms in anatase. 


Rutile-—-Schréder’s measurements of = in this crystal are confined 


to three wave-lengths, viz., 5461, 5790 and 6234A. These three readings 
exhibit such large casual errors that they cannot be used for a computation 
of thermo-optic behaviour. 


In general, neither the dispersion nor the thermo-optic behaviour of 
the three forms is known with sufficient accuracy to enable one to calculate 
the dispersion frequencies and their temperature shifts in a systematic 
manner. The data for anatase and brookite show, however, that the theory 
of thermo-optic behaviour for anisotropic crystals. which postulates a 
transfer of transition probability, as well as a temperature shift of the dis- 
persion frequencies, can explain these cases. The transfer terms are found 
to be very large and contribute substantially towards the temperature 
variation of birefringence. On the cther hand, the values of X are com- 
paratively small. 


4. THE REFRACTIVE INDICES OF RUTILE 


An absolute calculation of the birefringence of rutile and anatase has 
been made by Hylleraas (1926), based on Born’s method (1918). We shall 








— = c@ \y =) 
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perform this calculation in the case of rutile by the much simpler method 
of Bragg (1924), to see whether the very large refraciive indices and 
birefringence are explained in terms of the influence of crystal structure on 
polarisability. The rutile structure is tetragonal with the basis Ti (000; $44) 
and O(+ ww0; $+ w, 3-- ow, 1), with w=0-31, e=4-58A and 
¢ =3-73A. The structure can be pictured this way. The titanium atoms 
are at the corners and the body-centre of the tetragonal unit cell. To each 
Ti at the corners are attached two oxygen atoms at eaual distances, along 
the direction of one of the basal diagonals. To the Ti at the body centre 
are attached two oxvgen atoms at the same distances but in the direction 
of the other basal diagonal. 


In our calculations. it is more convenient to choose for the co-ordinate 
the two basal diagonals as the XY axes and the c or tetrad axis ss Z. 
The co-ordinates of the neighbours of an oxygen atom within a sphere of 
5A are easily calculated. Their description is given in Table VI, with 
reference to the chosen atom as origin. It will be seen that there are 42 
oxygen and 13 titanium atoms within the sphere. 

















TABLE VI 
Neighbours of an Oxygen Atom in Rutile 
Kind and No. . Distance from 
of atoms Co-ordinates origin (A) 
| ' 
| 2 Ti | wh oh 4 1-92 
}1 Ti os w w 0 2-07 
Lio eal 1-27 w 0 2-46 
18 0 se w dew 4 2-78 
120 ail 0 0 1 2-95 
1/20 ve) Qted 0 3-46 
| 2 Ti ae | w w-l 0 3-56 
2 Ti = w w 1 3-61 
12:0 a 1-2w 0 l 3-87 
10 cal 2w 0 0 4-00 
40 w| Qu} 3 1 4:55 
40 ial 4 } 0 4:58 | 
4 Ti = w wl 1 4-62 
8 O ool ze} l-w 4 4°87 
20 Pie 2w 0 1 4-97 | 
8 O w Sew 315 5-01 


From the mean index of the ordinary and extraordinary rays, one can 
calculate the molar refractivitv of O-- in rutile as 6-2 assuming that of Ti** 
to be 0-6 (cf. Fajans and Joos, 1924). Since the refractivity of the titanium 
ion is very small, one may assume that it is polarised as if in an isotropic 
medium, The polarisability tensor of the oxygen atom can be evaluated 
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by finding the field acting on it in different directions caused by the 42 
oxygen and 13 titanium neighbours. The calculations are considerably 
simplified because the cross-products xy, yz, zx vanish by virtue of the sym- 
metry of the crystal. Thus there is no component of the internal field at 
right angles to the electric vector. 


When the electric vector is parallel to the optic axis, the calculations 
are very straightforward. When it is perpendicular to the optic axis (say 
along Ox or Oy) the polarisabilities of the oxygen atoms in the layers 0, + ¢, 
+2c, etc., will be different from those of the atoms in the layers 


+ 4, + 5, -.--. But these values of polarisability are exactly interchanged 


when the electric vector rotates through 90°, by virtue of the tetragonal 
symmetry. Thus the refractive index is isotropic for directions at right 
angles to the optic axis. It is sufficient to calculate the polarisability of 
either type of oxygen atom for both the x and y directions of the electric 
vector. The mean of these values will be equivalent to the mean polaris- 
ability of both types of atoms for either direction of the electric vector. Thus 
the calculations become considerably simplified. 


Bearing these points in mind we can write for the refractive indices n, 
and n, (n, and n,) 
nii— l : p . 


where M is the molecular weight and p the density (4-22) of rutile while 


Ry and R,, are the molar refractivities. The anisotropy factor c, is given 
by Bragg’s theory as 


ont + Gul FCEE) eFC) me] 


1 





where €*= z* or (x*+ v?)/2 according as k is 3 or 1 and N is the Avogadro 
number. On calculating the summation in the right-hand side of (6), it 
turns out that c;= 1-030 and ¢,== 0-965. 


Substituting these values in Eq. (5) the two refractive indices for the 
sodium line are found to be as shown below: 
Calculated Observed 
2°86 2-89 


n 2 
Me 2:61 2-62 


The calculations show that these high indices and birefringence are 
caused by the large refractivity of the oxygen atom in rutile (nearly double 
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the value found in other crystals), as well as by a strongly anisotropic 
polarisability tensor (as in calcite) for instance 


Similar calculations made for anatase taking into account the influence 
of all neighbours within a sphere of radius 6 A gave the correct sign for the 
birefringence, i.e., 7, <M, but the magnitude of the calculated birefringence 
was twice as large as the observed value. The author believes that more 
neighbours will have to be considered if the birefringence is to be determined 
accurately both in the case of anatase as well as of rutile, and the excellent 
agreement found with rutile is fortuitous. 


T express here my thanks to Professor R. S. Krishnan for his kind 
interest in this work and to Dr. G. N. Ramachandran for suggestions and 
advice. 

SUMMARY 


Dispersion formule valid in the visible region are offered for rutile, 
anatase and brookite. These incorporate a characteristic frequency at 
2870 A in all the cases. The value of dn/dt in the three crystals are explained 
by large transfers of transition probability and moderate shifts towards 
larger wave-lengths of this frequency as the temperature is raised. In addi- 
tion, a calculation is made of the refractive indices and birefringence of 
rutile on the basis of Bragg’s theory and the results are in good accord with 
the observed data. 
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1. INTRODUCTION 


THE mechanical birefringence in liquids (Maxwell Effect) was extensively 
studied by VGrlander and Walter (1925). This led Raman and Krishnan 
(1928) to formulate a molecular theory for this phenomenon, in analogy 
with their theory explaining the electrical and magnetic birefringence of 
liquids. A quantitative comparison of the theoretical values with the experi- 
mental results was attempted by Paldhikar (1936) and Rao (1937) who 
showed that the calculated birefringence in the few cases verified by them was 
of the right order of magnitude. Since these papers were published, 
Buchheim, Stuart and Menz (1939) experimentally studied the Maxwell 
Effect in a few more simple organic liquids, thus providing a wider field for 
the verification of the theoretical ideas on the subject. An attempt is made 


now to critically examine these results in the light of Raman and Krishnan’s 
theory. 


2. RAMAN AND KRISHNAN’S THEORY 


Raman and Krishnan have shown from theoretical considerations, that 
the mechanical birefringence is a universal property which liquids may be 
expected to exhibit provided they are both geometrically and _ optically 
anisotropic. The Maxwell constant V is shown to be equal to 

(n®? — 1) (n? +2) 2 {(a — b)}(A — B) 
~WOnvKT ~~ (Say x (JAY? (1) 
where n= refractive index of the liquid, 
v = number of molecules per unit volume, 
K = Boltzmann Constant, 
T = absolute temperature. 


a, b, c are the principal geometrical dimensions of the molecules and 
A, B. C are their principal optical! polarizabilities. 


For a molecule possessing an axis of symmetry with @ = 4, equation (1) 
may be written thus; 
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(n?— 1)(m?+ 2) 2(c—a) ,, 5 
= -+. — - —— ° ‘Os & 
"= 10 nvKT Ja+c ¥ (2) 
where 6 is the optical anisotropy. 
3. DATA AND CALCULATIONS 


Equation (2) may more conveniently be put down thus: 

_ (n?— 1) (n?+ 2) M 2(c — a) 

= n “JOTRd 2atec 

The several constants given in Table I except a, c and 5 for different 

liquids are obtained from International Critical Tables. The geometrical 

dimensions of the molecules, collected in Table I, are generally accepted 

values, arrived at from X-ray data and other evidence. Even though slight 

geometrical asymmetry in the plane is obvious in the case of some of the 

benzene ring molecules in Table I, it is not taken account of in putting down 

the dimensions, as this difference may be ignored in the first approximation. 

The value of & for the liquids 1-* is taken from the existing literature 

wherever possible. For two liquids’® it is calculated from a knowledge 
of 8, and p. 


J8 (3) 


5 for mesitylene is not available from experiment and is therefore 
estimated by arguing that the structural change is progressive from benzene 
to tcluene to m-Xylene to mesitylene. 5 for toluene is 20-4 x 10-3, for 
m-Xylene it is 18-8 x 10-* and hence it is reasonable to extrapolate and 
obtain 6 for mesitylene as 17-2 x 10-* (i.e., the difference between the optical 
anisoptroy of m-xylene and that of mesitylene is assumed to be the same 
as that between that cf tcluene and of m-xylene). Evidence to this end is 
also forthcoming from the fact that as we go from benzene tc toluene, the 
symmetrical structure of the ring is disturbed: but the symmetry is pro- 
gressively reclaimed as we go from toluene to m-xylene and then to mesitylene. 
Thus mesitylene which is more symmetrical than m-xylene and toluene 
should have its 8 in the vicinity of that of benzene. 


8 for o-xylene is also not available from experiment and is estimated 
by arguing that the adiabatic compressibilities of the three xylenes are 
nearly the same (about 63 x 10-*) and therefore, the isothermal com- 
pressibilities for the three xylenes should also be nearly the same. 8, for 
p-xylene is 92 x 10-® and 8, for m-xvlene is 91 x 10-®. Hence it is reason- 
able to assume that 8, for o-xylene is 90 x 10°. The depolarization 
factor » for o-xvlene is 40%. Using these values of B; and », % is calcu- 
lated in the usual manner. * 

The figures thus collected and the results of calculation, using formula 
(3), are shown in Table I. Experimeritally observed values of the. Maxwell 
Constant are compared with the corresponding calculated results, 
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Liquid 





P-Xylene 
Mesitylene 
M-Xylene 
Chlorobenzo! 
O-Xylene 
Toluol 
Benzene 
Heptyl alcohol 
Carbon tetra- 
chloride 
Cyclo-hexane 





TABLE I. Maxwell Effect in Liquids 
| | Maxwell | Maxwell | 
.. Opaies al | 
| Motecutas| |netrcrive | Length of constant _ constant | 
weight | Density | | index | | Temp.) Anise | aces. A.U. | os experi- | Refer- 
| | tropy acid mental v7 lence for 
| eer 5 
j | | | 
um |e |» Ite. | xto? | oe | ¢ | elo? | xt 
i \ i 
mp pe _— — 
106-1 | 0-861 | 1-496 | 20 | 38-5 7 5 164-0 275-2 (7) 
120-1 | 0-863 | 1-497 > | S88 9 3 122-8 196-1 oes 
106-1 | 0-865 | 1-497 vs Is-S 7 : 115-5 193-2 (1) 
112-5 | 1-107 | 1-526 » | 30-6 7 ; 130-4 | 186-1 (4) 
106-1 | 0-879 | 1-506 «| SES 7 $ 86-8 | 183-7 a 
92-06 | 0-866 | 1-496 » | 20-4 7 } 104-0 155-6 (2) 
78-05 | 6-878 | 1-501 >» | 180 7 ; 2+7 96-1 (3) 
116-1 | 8-817 | 1-425 ‘i 2-2 4-6 | 11-34 52-37 58-43 (5) 
153-8 1-595 | 1-46! a. 1-6 a a <8+7 (6) 
84-09 | 0-779 | 1-427 1+8 8 5 15-7 | <43 (8) 
} | 
4. Discussion OF RESULTS 
The following observations may be made in respect of results in 


Table I: 
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In all the liquids, there is a fair agreement between the theoretical and 
experimental values—in fact, the agreement is much better than that shown 
in the earlier work by Paldhikar (1936) and Rao (1937). 


In all the cases, except cyclohexane, where the very low experimental 
value might not have been measured with reasonable accuracy, the theoretical 
value is less than the experimental value. 

The ratio of the theoretical value to the experimental one in most cases 
is, surprisingly enough, almost the same, roughly 4. 

The larger experimental value cannot be attributed to turbulence cr rise 


of temperature. 


to be correct. 
In conclusion, the author desires to express his grateful thanks to 


Nevertheless, the main ideas underlying the theory appear 


Dr. S. Bhagavantam, Director, Physical Laboratories, Osmania University, 
Hyderabad, Deccan, for kindly suggesting the problem and for his keen 
interest and encouraging guidance in the course of the preparation of this 
paper. 
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1. INTRODUCTION 


A POLYCRYSTALLINE aggregate will be elastically isotropic when the crystals 
that compose it are distributed uniformly over all directions, and in such a 
case, it will exhibit only two elastic constants which are independent of 
direction. Bridgman (1931 and 1932) determined the compressibilities of a 
number of cubic compounds by a static method. He worked with poly- 
crystalline specimens obtained from dry powders after compression to 
about 310,000 Ib. /in.2 in a mould. Lawson (1940) also determined the com- 
pressibility, among other properties, of ammonium chloride using the 
composite piezoelectric oscillator method. He obtained the polycrystalline 
specimens from a block of material prepared by compressing dry granulated 
ammonium chloride to 50,000 Ib./in.? in a steel cylinder at room temperature. 


The ultrasonic wedge method (1944) has been successfully used by us 
to study the elastic properties of single crystals during the past several years. 
It is now proposed to extend this method to a study of the elastic properties 
of polycrystalline aggregates. The present paper embodies the results of 
a preliminary investigation of some cubic compounds. There is a specia] 
advantage in this method, because it can be used for investigations at higher 
and lower temperatures without much difficulty and this will be attempted 
in our future work. 


2. MATHEMATICAL THEORY 


The calculation of the properties of polycrystalline aggregates is based 
on the idea that the properties of the aggregates are average values of the 
properties of the individual crystals. The average value is obtained by 
summing up the values of the property of the crystals according to the 
frequency distribution of their orientations, and dividing this sum by the 
number of orientations employed. In the case of random orientation, the 
crystallographic axes are distributed with uniform frequency in all direc- 
tions of space. Details of this process are given in a recent publication of 
Boas (1949), 
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In this way, the average value of C,,’ can be obtained from the relation 
us 27 


J C,,’ sin 0 dé db 


(Cuy')ay = 0-9" 


T 2 


f f sin 0 dé db 


0 
where C.. = (2C,,+ Cy) ot (C,i- Cio— 2 Cua) (a,54+ Ogg G9") in the case 
of cubic system. 3, ay3 and a3, are the direction cosines and are given by 
a43= Sin 8 cos d, ag3= sin 8 sind and az,= cos @. 
Working out the integration, the expression for (C,,’),4y comes out as 
; 3 _.2 
(Ci), dav = s Ci. 4 5 (Cie +2 C44). 


The average value for C,,’ can be obtained by performing the integra- 
tion 


27 


db ( db (Ca' sin 6 dé, 
(Caav = f om nf — 
vu 0 0 


where Cqq’ = (049? a43?-+ Gag” 937+ agg” agg”) (Cyy— Cyga— 2 Cag) + Cay for 
cubic system. 


443, %3 and agg are the same as above and 
d42= — cos 8 cos ¢ sin s — sin ¢ cos #; 
Aoo= COS ¢ CoS # — sin d sin % cos 8 and a;.= sin 4 sin wf. 


Introducing these values and working out the integration, the expression 
for the average value for cubic crystals becomes 
, 3 1 
(Caa)av = 5 Cag t+ 5 (Cu — Cy). 


These expressions have been derived also by Mason and McSkimin 
(1948 and 1949) in order to calculate the scattering factors for metals. 


It may be noted that (Cy,’)ay — Sihe = = which is the 


bulk modulus of the polycrystalline aggregate. This is not a third inde- 
pendent elastic constant. Further, if the crystal is such that the Cauchy 
relation C,,= C,, holds good, the above equations become dependent on 
one another and reduce to 


, , 3 
(Ci; Jav =3 (Cy, dav = 5 (Ci; +2 Cy). 
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From this, it follows that the bulk modulus K is given by the relation 
5 , 5 ’ 

K = 9 (Cy )av = 3 (Caa’av- 


3. EXPERIMENTAL TECHNIQUE 


The experimental set-up used is essentially the same as that employed 
in the wedge method. Some of the alkali halides belonging to the cubic 
system, the elastic constants of which have been determined in the crystalline 
state, are chosen for the present investigation. Pure and dry substances 
are obtained, finely powdered and passed through a 200-mesh sieve. These 
powders are compressed in a steel mould between two high tensile steel 
cylinders of 1” diameter. All the polycrystalline specimens in the present 
investigation are pressed at about 60,000 Ib./in.2, The compressed cakes are 
then ground to uniform thickness on very fine emery cloth. Thus a uni- 
formity of 1 to 2 in 100 is obtained. That the microcrystals composing the 
polycrystalline aggregate are oriented in a perfectly random manner as a 
result of such a compression has been tested by an X-ray study of one or two 
specimen tabloids. 


The determination of the densities of these compressed cakes is a very 
important factor since the calculation of (C,,’),y and (C44’),y involves a 
knowledge of the densities besides the appropriate characteristic frequencies. 
In all cases, the density is determined directly by evaluating the volunie of 
the circular tabloid and determining its weight in air. In this way, densities 
of the various tabloids have been determined to an accuracy of 1 to 2 per cent. 


It is shown by Walker (1923) that when a powder is pressed into a 
cylinder by a plunger, the density of the powder is not evenly distributed 
but is greatest immediately beneath the plunger, gradually falling off as the 
distance from the plunger increases. Attempts to detect such large variation 
in density have yielded negative results in the case of the present tabloids, 
because the maximum thickness does not exceed 3-5 mm. and in some cases, 
the thickness is only about 2 mm. 


In the present investigation, while determining the characteristic fre- 
quencies of the tabloids, special care is taken to see that the tabloids do not 
come into direct contact with the liquid in the cell. Otherwise, the liquid 
gets into the pores of the specimens and vitiates the experimental results. 
For this purpose, the design of the crystal holder is slightly modified. The 
lower annular brass electrode on which the specimen under investigation 
rests, is now replaced by a brass cup-like vessel. The central hole of the 
latter is covered by a thin aluminium foil. The tabloid is placed on the 
aluminium foil with a drep of transformer oil in between for securing good 
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acoustical contact. The quartz wedge is kept on the top of the tabloid with 
another drop of transformer oil in between as usual. This oil, being dense, 
does not enter the pores of the tabloid as has been revealed by the constancy 
in the mass of the specimen just before and after the experiment, the weight 
being determined after wiping out the oil layer on the surface. But for this 
minor modification in the design of the crystal holder, the rest of the experi- 
mental set-up is just the same as that used for single crystals. The charac- 
teristic frequencies of the compressed specimens are determined in the usual 
way in the frequency range 2 to 10 Mc/Sec. The frequencies are measured 
using a precision heterodyne wavemeter capable of reading correct to 
within 1 in 50,000. 


A feature to be mentioned here is that the power of the oscillator 
exciting the polycrystalline aggregate should neither be too great nor too low. 
If it is too great, a continuous diffraction pattern consisting of forced 
oscillations is observed instead of discrete maxima; if it is too low, no 
maxima will be observed due to the excessive absorption of the ultrasonic 
beam by the polycrystalline aggregate. Hence these two factors have to 
be reconciled and an optimum power just sufficient to give the discrete 
maxima should be used. 


4. EXPERIMENTAL DATA AND RESULTS 


The densities of the polycrystalline specimens pressed to about 
60,000 Ib. in.*, determined by the method outlined in the previous section, 
are given in Table I along with the densities of the corresponding single 
crystals. 

TABLE | 
Densities of compressed tabloids 


No. | Substance Prabloid g.iem.” | Perystal g./em.* 
] | KCl 1-96 1-984 
2 NaCl 2-12 2-168 
3 NH,Cl 1-50 1-527 
4 | KBr 2-67 2+ 756 
5 | KI 2-95 3-127 
t 


| 


The average values (C,,"),y and (Cq4’),y of the elastic constants calcu- 
lated from the characteristic frequencies and the observed densities of the 
tabloids (column 3 in Table I), are given in units of 10% dynes/cm.? in 
Table If. The values of (C,,')sy and (Cq,’),y calculated from the known 
single crystal constants using the expressions given in the second section, 
are compared with the measured results. 
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TABLE Il 
| | | | ae 
a is | Value 
ae _ | Fundamental Effective 
No.| Substance —— | frequency = elastic |— - gowneetnte 
| | in Mc./Sec. | constant | Measured | Calculated 
| j | i) eee” 
-\2EO eo} = 9099 0-866 Vg (Car Jav 3-08 3-18 
do... | do 0-501 t (C44’)av 1-03 1-05 
2 | NaCl ..' 3-24 0-705 Z OSes ten 4-42 4-51 
do weil do 0-407 t (Cee hav 1-47 1-50 
3 | NH,CI ..| 2-89 0-785 / CyJav | 3-09 3°17 
| do | do 0-448 t ( 44 av | 1-01 1-04 
4 | KBr | 3-00 0-5¢3 Z Cale | oe 2-73 
| do do 0-290 t (Cale | = OO 0-92 
5 | KI «| 188 0-644 d (C0 daw | 164 | Ql 
do. do 0-378 t Cake | On |lCUe 
} ; 








* Here, / and ¢ denote the longitudinal and transverse modes respectively. 
The values of the bulk modulus K calculated from the above experi- 


mental results using the relation K = atti hee are compared in Table III 


with those measured by direct experiment, in units of 10" dynes/cm.* 








TABLE III 
Sabstance | (Cu‘Jav | (Caa’Dav 0. Geos dame 
| 
KCl ..| 3-08 1-03 bm} oT 
NaCl re 4-42 1-47 2-46 2-40 
NH,Cl ..| 3-09 =| = 1-01 1-72 1-79 
KBr we] 2048 0-81 1+35 1-50 
KI es 1-64 0°57 0-91 1.17 





The experimental data concerning ammonium chloride and magnesium 
oxide available in the literature are given in Table IV. 








TABLE IV 
7 T 
Substance Porystal | tabloid K Observer 
: 7 
NH,Cl | 16527 és (097 Sundara Rao and Balakrishnan 
| 1.79 | Bridgman 
1 +520 1-77 Lawson 
1-512 1-66 do 
MgO ee 3°57 oe 16-6 Bridgman 


do 
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5. DISCUSSION OF RESULTS 

Table I shows that the densities of the polycrystalline specimens 
approach the corresponding single crystal densities to within about 2 per 
cent. only in the case of the first three substances, viz., potassium chloride, 
sodium chloride and ammonium chloride. The tabloid densities for the 
bromide and iodide of potassium are lower than the corresponding crystal 
densities by more than 3 per cent. and 6 per cent. respectively. These results 
show conclusively that all substances cannot be compressed to the same 
extent by the application of the same pressure. Similar conclusions were 
reached by Hagen (1919) who showed that it is not possible to compress 
quartz and barium sulphate to cohesive tablets while the alkali and. other 
halides give very cohesive products. Further, the fact that the density of 
the pressed specimen of magnesium oxide with which Bridgman (1932) 
worked was only 3-01 in spite of the large pressure employed by him while 
the crystal density is 3-57, also shows that this substance is even worse than 
potassium iodide from this point of view. 

Table Il shows that in the case of the first three substances, whose 
densities lie within 2 per cent. of the crystal density, the measured average 
values of C,,’ and C,,’ are also lower than those calculated from the known 
single crystal constants by only about 2 to 3 per cent. In such cases, we 
may extrapolate the density and obtain the ideal elastic constant if required. 
On the other hand, in the case of potassium bromide and iodide where the 
densities of the specimens are lower than the crystal densities by more than 
3 per cent. the deviations in the average values of elastic constants are 
disproportionately larger. In other words, the compressibility of a pressed 
specimen increases rapidly and out of proportion to the decrease in density, 
if we depart appreciably from the ideal densities. These results are in agree- 
ment with the observations of Bridgman in the case of magnesium oxide 
(Table 1V) where a 18 per cent. fall of density is followed by a 40 per cent. 
fall in the bulk modulus. In our work, a 6 per cent. fall of density in the 
case of potassium iodide is followed by a 25 per cent. fall in the bulk modulus. 
On the other hand, Bridgman (1931) has also determined the compressibilities 
of a number of cubic substances like NH,Cl, NaCl, KCl, KBr, etc., pressed 
at about 310,000 Ib./in.2 and his values of compressibility agree well with 
those for single crystals available in the literature. But, as Lawson remarked, 
the densities of these specimens had not been reported by Bridgman and 
we have to presume that they are very close to the respective single crystal 
values. 


In the case of ammonium chloride, Lawson however observed, as the 
figures in Table IV reveal, that the isothermal compressibility of a specimen 
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whose density is 1-512 gm./cm.* is 6-7 per cent. greater than that of a speci- 
men whose density is 1-520 gm./cm.? It may be seen from the data that 
Lawson’s value for the bulk modulus for a specimen of density 1-520 gm./cm.* 
is slightly less than the value obtained by Bridgman. The latter author, 
however, did not give the density of the specimens with which he worked. 
Sundara Rao and Balakrishnan (1948) determined the elastic constants of 
single crystals of ammonium chloride. Their values lead to a bulk modulus 
of 1-77 x 1021 dynes per cm.? which is in good agreement with the isothermal 
value of Lawson. However, Lawson’s experimental result, that the lowering 
of density by 1 per cent. causes a fall of 6-7 per cent. in the bulk modulus 
does not fit in with the results of the present investigation; nor is it in agree- 
ment with the experimental results of Bridgman in the case of magnesium 
oxide. 

The bulk moduli calculated from the values of (C,,’),y are compared 
in Table I[L with the standard values available in the literature. For the first 
three substances in which the densities approach the ideal values to within 
about 2 per cent., the moduli also approach the single crystal values to within 
about the same percentage. 


6. SUMMARY 


The ultrasonic wedge method has been used for a determination of the 
average values (Cyy")ay and (C4,4’),y in the case of certain polycrystalline 
aggregates of cubic compounds prepared by compressing finely divided 
powder under high’ pressures. The measured values are compared with the 
values calculated from the known single crystal constants. It is found that 
certain substances can be compressed till their densities closely reach the 
single crystal densities and in those substances the average elastic properties 
also closely reach the corresponding ideal values. In other cases, the com- 
pressibility of a pressed specimen depends to a marked degree upon its 
density. 
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INTRODUCTION 


SEVERAL methods have been described in literature to measure the variation 
of boundary tension with time. Adam! has employed the maximum 
bubble pressure method to study the variation of surface tension with time, 
for a period of ten minutes, but concludes that the method is not suitable. 
Venkatachala* and Krishnappa® arrive at the same conclusion after employ- 
ing this method to study the ageing of surfaces of aqueous solutions of wetting 
agents. The sessile bubble method has been used by Adam and Shute! for 
studying the variation of surface tension with time, of certain paraffin chain 
salts. The method seems to be free from contact angle effects but special 
precautions have to be taken against accidental contamination. Tartar, 
et al.* have made a comparative study of the capillary rise and the sessile 
bubble methods for the determination of the variation of surface tension. 
They prefer the latter method and claim that it gives reproducible results. 


The trough method suggested by Doss,* and McBain and Wilson,’ is 
quite suitable for a study of the variation of surface tension with time but 
this method is cumbersome for the study of the variation of interfacial 
tension with time. The pendant drop method described by Andreas, 
Hauser and Tucker, appears to be satisfactory for the latter purpose, but 
the apparatus is very complicated. 


In the present method, the advantages of the pendant drop method 
are combined with those of the drop weight method. The method consists 
in forming a small drop of the solution at the end of a vertical capillary 
placed inside a thermostat. The drop is allowed to age, and the time 
required for the drop to fall, is noted. The boundary tension is calculated 
by applying Harkin’s* corrections for the volume of the ideal drop. 


EXPERIMENTAL 
Reagents empolyed 


(i) Benzene.—The benzene (Kahlbaum’s pure grade) was distilled in 
an all-glass apparatus and the middle fraction was used. 


* Part of the Thesis submitted by R. S. Subrahmanya for the M.Sc. Degree of the 
Mysore University. 
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(ii) Congorubin.—Congorubin was purified by the sodium acetate 
method of Robinson.“ It was washed with alcohol and dried in vacuum. 

(iii) Potassium iodide, acetic acid and sodium chloride-—A. R. quality 
reagents were used. 


Apparatus 


The stalagmometer used for the measurement of boundary tension 
consisted of a calibrated micro burette which had a suitable capillary tube 
fused on to its delivery end. The tip of the capillary was carefully ground, 
using fine corborundum powder,’ 1° and finally wet talc. As pointed out 
by Andreas, et al.,® microscopic perfection was not necessary since the liquid 
tended to bridge slight irregularities. The capillary tube passed through a 
rubber stopper fixed to a small bottle containing benzene, so that the tip was 
immersed about 0-5” below the surface of benzene. The stopper carried 
a sensitive thermometer and a small capillary tube, open to the atmosphere. 
The entire apparatus was placed in an air-thermostat at 30°C. (+ 0-05°). 
The radius of the capillary was 0-3402 cm. 


Determination of Interfacial Tension between Water and Benzene 


The stalagmometer was cleaned thoroughly’ with warm chromic acid 
and finally washed with conductivity water. It was then filled with con- 
ductivity water saturated with benzene. Ne grease was employed to 
lubricate the stopcock. The bottle was filled with benzene saturated with 
water, and kept in the thermostat for five to six hours. The stopcock was 
gently opened and a small drop of water was allowed to form in the benzene. 
The drop was sucked back so that the liquid meniscus was at a fixed mark 
on the capillary tube. The water level in the micro burette was then noted 
by means of a travelling microscope. The stopcock was next opened and 
a fresh drop of water in benzene, was allowed to form slowly. The time 
of formation of each drop was about five minutes. After the drop fell. off 
from the tip, the water was once again sucked back to the mark, and the 
stopcock closed. The microburette was then read and the volume of the 
drop determined. Within the limits of experimental error, it was noticed 
that the volume of the drop of water was independent of the time of forma- 
tion of the drop. Evaporation of water from the microburette was mini- 
mised by attaching to its open end a fine capillary. The reproducibility 
in the formation of drops can be gauged from the following values obtained 
in an experiment :—0-3903 c.c., 0-3909 c.c., 0-3903 c.c., 0-3887 c.c.—mean 
value 0-:3901 c.c. The interfacial tension y was calculated from the following 
formula :— 
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Vv (D, 2a D,)g 
ie 2R6-” 


where, R = radius of tip of capillary, 
V =volume of drop, 
D.= density of water, 
D,= density of benzene, 


g = acceleration due to gravity. 


R/V" being known the value of ¢ was ascertained from data furnished 
by the Harkins and Brown.® Using the above formula, the interfacial tension 
was calculated. The mean value for the benzene-water interface was 34-29 
dynes per cm. which agrees with the standzard value of 34-34 at 30° within 
the limits of experimental error (about 0-3%). 


Determination of the Variation with Time of the Interfacial Tension between 
Benzene and Aqueous Solutions of Rubin 


As pointed out already, the chief difficulty in the determination of the 
variation of interfacial tension with time, was the maintenance of the inter- 
face in an undisturbed condition during the entire period, over which the 
measurement extended. In the present investigation, this difficulty was 


overcome by employing the principle of the pendant drop method. The 
following procedure was adopted :— 


After cleaning, the stalagmometer was filled with the rubin solution 
saturated with benzene, while the bottle was filled with benzene saturated 
with the rubin solution. After the attainment of temperature equilibrium, 
a drop of known volume was produced in four to five seconds at the end of 
the capillary tube. The drop was allowed to hang on at the tip of the 
capillary when it was noticed that the semi-spherical shape of the drop 
gradually changed to that of an elongated sphere until finally the drop fell 
into the benzene. The interval between the formation and the falling of the 
drop was noted. The aqueous solution was sucked back to the mark and 
the final volume of the drop determined. The interfacial tension at the 
time of falling of the drop was determined by Harkin’s method, as already 
described. It was noticed that for a given concentration of rubin, the time 
of falling of a drop was a function of the volume of the drop. The intervals 
for the falling of drops having various volumes were noted and the corres- 
ponding values of interfacial tensions calculated. Fig. | represents the 
variation of interfacial tension with time for different concentrations of rubin 
solution. The rubin solution (0:005%) was then mixed with a few drops 
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Fic. 1. Effect of concentration on the variation of interfacial tension, with time, of 
Congorubin solution at benzene interface. 


Concentration of congorubin pertaining to curve 1= 0-00005%, curve 2 =0-0005%, 
curve 3 = 0:005%, and curve 4 = 0-05%. 
of dilute acetic acid to lower the pH to 3-8 and fresh measurements made 
with the solution. From the results obtained (Fig. 2) it is clear that a lower- 
ing of pH greatly diminishes the interfacial tension. ' 


Effect of Salts on the Variation of Interfacial Tension with Time 


Evans!% and Guest!* have shown that salts of the type of sodium or 
potassium chloride increase the interfacial tension between aqueous solutions 
and an inert organic liquid, while potassium iodide decreases the tension. 
It was of interest to examine if sodium chloride and potassium iodide have 
the same effect on the interfacial tension between congorubin solution and 
‘benzene and how the interfacial tension varies with time. 0-005% rubin 
solution was employed. The salt concentration was 0-1 molar. The 
results obtained are indicated in Fig. 3. It will be noticed that both the 
salts considerably lower the interfacial tension, the effect being greater with 
the iodide. With time, the interfacial tension decreases in both the cases. 


Measurement of Surface Tension 


It was found necessary to prevent evaporation of water from the 
pendant drop, by having in the bottle some of the rubin solution, 1 to 2mm, 
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Fic. 2. Effect of pH on the variation of interfacial tension, with time, of Congorubin 
solution (0-005%) at benzene interface. 
below the drop. The variation of surface tension with time for the rubin 
solution was also determined by the trough method. The results obtained 
by the two methods are given in Fig. 4. The agreement between the values 
obtained py the two methods is very satisfactory. 


DISCUSSION 
Variation of Surface Tension with Time 


Investigations by the present method, on the variation of surface tension 
with time, using 0:2% solution of congorubin in water, confirm the results 
obtained by the trough method. The present method can be employed for 
a study of surface tension at different temperatures more easily than the 
trough method. Another advantage of the method is that only small quanti- 
ties of solutions are required. The slow variation of surface tension with 
time shows that the accumulation of the rubin molecules at the surface of 
the solution is of the activated type and is a slow process (cf., Doss.1> McBain, 
et al.**), 


Variation of Interfacial Tension with Time of Aqueous Rubin Solutions at 
Benzene Interface 


The effect of concentration of rubin on the variation of interfacial ten- 
sion with time is quite marked. At very low concentrations (0-00005%) of 

















Method for Study of Variation of Boundary Tension with Time 141 


a 
/ 


30F 








Interfacial Tension in Dynes per cm. 
& 











S 1 | 


a a a ve 


Time in Minutes 





( a 2 
, 5 


Fic. 3. Effect of Salts on the variation of interfocial tension, with time, of Congorubin 
solution at benzene interface. 


Composition of aqueous solution :— 
Curve 1 0-005% dye 
oo 2 . +0:1N sodium chloride. 
- 3 i +0:°1N potassium iodide. 
» 4 O-1 N potassium iodide. 


(Note.—As with potassium iodide solution, there was no change in interfacial tension 
with time, in the case of 0-1.N sodium chloride). 
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rubin, the interfacial tension will be practically the same as that between 
water and benzene. There is scarcely any change with time. But at con- 
centrations of the order of 0:0005% and above there is not only a fall in the 
interfacial tension but also a variation in the tension with time. This change 
is rapid in the beginning but gradually falls off. At the air-solution interface 
on the other hand, there is no variation in surface tension during the first 
20 minutes (Fig. 4). The tension then begins to diminish and we finally 
get a practically constant value which is characteristic of a condensed film. 
Lower concentrations (< 0:2%) of rubin cannot be employed for the study 
of surface tension owing to very slow accumulation. In studying interfacial 
tension, however, solutions of far lower concentrations (0-05%) can be 
used. A fall in interfacial tension with time is noticed in such solutions. 
The film in this case is of the gaseous type as can be gathered from the nature 
of the curve, while the corresponding film for the air-solution interface is 
of the condensed type. The difference can be explained thus :— 


Congorubin has four hydrophilic groups as shown by its constitution. — 


—— 
CoH — N= N — [2] ~ Cols (4) So,Na 


| 9) 
Ci, — N= N—pi—CGe | i 
As shown by Doss?” for benzopurpurin, the slow surface accumulation 

of rubin may be due to the extra activation energy required to eject two of 
the hydrophilic groups into the air. Thus the activated molecule will have 
the hydrophobic group and also the hydrophilic group pointing out towards 
air. This explains the formation of the condensed film at the air-water 
interface. But when air is replaced by benzene there will perhaps be no 
tendency for the two hydrophilic groups to come out. On the other hand 
the hydrophilic groups will all be pointing out towards water while the 
hydrophobic portion will get oriented towards benzene, as in the case of 
sodium oleate solutions.» Hence the only force that retards the rate of 
accumulation at the benzene interface is the electrical double layer formed 
at the interface due to accumulation of rubin molecules. Owing to the 
potential of the electrical double layer, the concentration of the capillary 
active ion, in the region just beneath the surface layer, will be very much 
suppressed and will be equal to™® 


Eygr 
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where c = bulk concentration of the solution, 














t 
' 
t 








Method for Study of Variation of Boundary Tension with Time 143 


e = charge on the electron, 


thy= potential due to the surface fully covered by the solute molecules, 
x = fraction of the total surface covered, 

k = Boltzmann constant, 

T = absolute temperature. 


Thus, accumulation at the surface occurs from a region of low concentra- 
tion. Hence the rate of accumulation will be proportional to 


eysoX 


e~ F- 

Examination of this equation shows that the rate of accumulation should 
be very high at the start (i.e., for small values of x). As the surface gets 
enriched with respect to the solute molecules (i.e., x increasing) the rate of 
accumulation goes onjdiminishing exponentially. Thus the rate of fall in 


interfacial tension will be high to begin with, but rapidly diminishes with 
time. 


Andreas® and others have determined the variation of surface tension 
with time, using sodium oleate solutions. They have also determined the 
variation of the interfacial tension between these solutions and mineral oil. 
Unlike congorubin, sodium oleate gives the gaseous type of films both at 
the air-water interface and at the water-mineral oil interface. Evidently, 
the accumulation of the sodium oleate is mainly governed by the electrical 
double layer at the two interfaces. Adam and Shute* have come to similar 


conclusions by working with cetyl pyridinium bromide and sodium dodecyl 
sulphate solutions. 


The rapid fall of interfacial tension (i.e., the high rate of accumulation 
of the rubin at the interface) in presence of electrolytes needs some explana- 
tion. Adam and Shute* have noticed similar effects with cetyl pyridinium 
bromide at the air-water interface. They state that the rapid accumulation 
is due to micelle formation. But the actual mechanism proposed is that the 
high negative charge on the micelles drives very rapidly the molecular species 
of the pyridinium bromide, to the interface. In the case of the congorubin 
solution-benzene interface the effect can be ascribed to the diminution in the 
value of 4%, (Zeta-potential) with an increase in the concentration of the 
electrolyte, in the rubin solutions. The rate of accumulation which is pro- 


portional to ce ve gets enhanced with the decrease in the value of y,. 
Hence the interfacial tension falls rapidly when the concentration of the 
electrolyte is increased. The effect of pH on the variation of interfacial 
fension with time can also be explained on the basis of the diminution of 
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the zeta potential. An increase in hydrogen-ion concentration reduces the 
zeta potential and thus increases the rate of accumulation (Fig. 2). 


SUMMARY 


1. A new method for the determination of variation of interfacial and 
surface tensions with time, has been described. In this method, the advantages 
of the pendant drop method are combined with those of the drop weight 
method, thus simplifying the experimental technique. 


2. Variation of interfacial and surface tensions with time has been 
investigated for congorubin solutions. The rubin solution yields at the 
air-liquid interface a solid film and at the solution-benzene interface a gaseous 
type of film. 


3. The effect of concentration of rubin on the variation of interfacial 
tension with time has been studied. At concentrations of rubin higher than 
0-0005% the fall in the interfacial tension will be large. 


4. Salts like sodium chloride and potassium iodide lower considerably 
the interfacial tension between rubin solution and benzene. Lowering of 
pH also brings about a fall in interfacial tension. 


5. The slow variation of interfacial tension has been accounted for, 
on the basis of changes in zeta-potential. 
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1. INTRODUCTION 


Tuis paper deals with the theory of some interesting phenomena which are 
observed in the transmitted X-ray beam when a Bragg reflection takes place 
in a crystal. Obviously, if the crystal is non-absorbing, then the whole of 
the incident energy must be transmitted in the absence of the Bragg reflection 
and when a reflection occurs, the energy which is reflected will manifest 
itself as a reduction in the intensity of the transmitted beam. Thus, the 
reflected and transmitted intensities would be complementary, so that if 
R and T are the reflection and transmission coefficients, then 


R+T=1. (1) 


If the crystal is absorbing, we have one more mechanism whereby energy 
may be abstracted from the incident beam. If A is the fraction of the 
incident intensity removed in this way, we have from a priori considerations, 
the following equation: 


R+T+A=1 or R+T=1—A=T, (say) (2) 


This means that T is always less than or equal to Ty, and near the Bragg 
reflection, when R has an appreciable value, there should be a reduction in 
the transmitted intensity. 


The foregoing deduction is, however, based on the idea that absorption 
and reflection are independent processes. However, if one supposes that 
the two are interrelated and that the absorption coefficient can itself be 
modified by the presence of an interference maximum, then conditions may 
be different. Thus, if the effective absorption coefficient is reduced for 
directions of incidence close to the Bragg angle, then the quantity A in 
Eqn. (2) would itself be reduced and one cannot conclude, without further 
calculations. whether the transmission would effectively decrease or increase. 
It is conceivable that the decrease in absorption coefficient would be such 
that the consequent increase in the transmited intensity may be larger than 
the decrease due to part of the energy being reflected. In such a case, one 
may expect a peak in the transmitted beam near the Bragg reflection rather 
than a dip. In fact, such a result has been observed by Campbell (1951) 
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using perfect calcite crystals. He found that the peak of the transmitted 
intensity may even be several hundred times as large as Ty. The name 
** anti-reflection * is suggested for this phenomenon. As is evident from 
the above discussion, such a phenomenon can be expected only with an 
absorbing crystal. If there is no absorption, Eqn. (1) would be rigorously 
valid and there would be always a diminution in T near a Bragg reflection. 
Also, the greater the value of A, i.e., the thicker the crystal, the more 
striking would be the intensity changes produced by a small alteration in 


the absorption coefficient, and consequently the more prominent would be 
the anti-reflection peak. 


These qualitative arguments are borne out by the detailed theory given 
below. It is based on the well-known dynamical theory of Ewald (1916). 
In particular, the development of the theory due to Laue (1949) forms a 
proper starting point for the discussion in this paper. Compact expressions 
have been derived for T and R, both in the so-called “* Laue” and “ Bragg ” 
cases, and these have been numerically worked out in the former case. They 
show clearly that, for increasing thicknesses of the crystal, the transmission 
curve slowly alters its form from a minimum to a sharp maximum near the 
Bragg reflection. 


2. PREVIOUS LITERATURE 


On the experimental side, the subject goes back to an interesting 
observation made by Rutherford and Andrade (1914) during their determina- 
tions of y-ray wavelengths using crystal diffraction. They observed in 
general a diminution in the intensity of the transmitted beam when a 
diffracted beam occurred. However, this diminution was not observed 
when a good perfect crystal was chosen for the experiment. In continuation 
of this experiment, Borrmann (194]) allowed a divergent beam of X-rays to 
fall on a crystal plate of quartz (0.2 mm. thick) and found lines in the trans- 
mitted beam where the transmission was appreciably reduced. However, 
when he chose a perfect crystal, as evidenced by the sharpness of the X-ray 
reflection, he found that the lines exhibited a ‘ Helldunkelstruktur’’, i.e., 
adjoining areas having an intensity larger and smaller than the general 
background. The occurrence of a transmitted heam of intensity larger 
than T, is particularly noteworthy. 


Lonsdale (1947) made similar studies with a number of crystals and 
confirmed Rutherford and Andrade’s observation. She found that the 
deficiency lines in the transmitted beam were imperceptible with a perfect 
crystal, while they became clearer with increasing degree of mosaicity. She 
did not make a critical study of the distribution of intensity of these lines, 
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Making use of the dynamical theory of X-ray diffraction, Laue (1949) 
has worked out the theory of the transmission cof X-ravs when interference 
takes place. In particular, he considered the so-called “ Laue case”’, i.e., 
when both the incident and diffracted beams emerge from the back surface 
of the crystal plate. Using his fermule he also calculated numerically the 
angular distribution of intensity of the reflected and transmitted beams in a 
few particular cases. These are similar to those observed by Borrmann 
(1941), although no quantitative comparison can be made. 


More recently, Campbell (1951 a, 6) studied the phenomenon quanti- 
tatively, using Geiger-Miiller counters. The experiments were performed 
with a monochromatised beam of Cu Ka having a small angle of divergence 
of the order cf 15 seconds of arc. Symmetrical *‘ Laue” reflections were 
obtained from the crystal under study and both the reflected and the trans- 
mitted intensities were measured simultaneously for different settings near 
the Bragg reflection by means of two Geiger counters. Sharp maxima were 
found both in the reflected and the transmitted beams with perfect crystals. 
With an imperfect crystal, the increase was imperceptible, but there was no 
decrease. Campbell has stated in his paper (1951 5) that the sharp maxima 
in the transmitted beam do not correspond to the curves given in Laue’s 
paper (1949) and that his results are therefore inexplicable on this theory. 
However, the conditions of his experiment are different from those assumed 
in the theory and when the proper conditions are incorporated in the theory, 
it leads to sharp maxima, as is shown below. 


3. NOTATION 


The notation follows mainly Laue’s symbology (1949), although a 
number of new symbols are used. Zachariasen (1945) has also given an 
account of the dynamical theory and his notation has been adopted in 
certain cases. The symbols within brackets are the equivalents in 
Zachariasen’s notation of those used here. 


A Wavelength of the X-rays in vacuum (Ag) 
k (= 1/A), length of the wave-vector in vacuum (Ky) 
h Index (subscript) denoting a particular X-ray reflection 

(This index may be omitted wher there is no confusion). (H) 
9 Angle of incidence (6) 
6, Geometrical Bragg angle (03) 


4, Bragg angle, allowing for refraction 


pe Linear absorption coefficient (in the absence of Bragg 
reflection) 


(4g) 
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D Thickness of the crysta! 


C Polarisation factor, = | for the normal component and = 
cos 2 4 for the parallel component 


Yo) Cosines of the angles made by the incident and reflected 
beams respectively with the normal to the surface 


Xo| Fourier components of order 0 and A in the triple Fourier 
xn) expansion of 47a, where a is the polarisability 


oat Real and imaginary parts of X,:X,= X,,+ iX;,; 


D, Electric vectors of the diffracted waves of order 0 and h 


Do } Magnitudes of Dy, D, 


D; 
Zz Unit vector, normal to the surface, directed inside the crystal 
mg Wave vectors, of order 0 and A inside the crystal 

h 
: } Subscripts to denote the two wave fields inside the crystal 
a 


Superscripts, a to denote the front, and d the back surface of 
d the crystal, e.g. 


D,'” represents the electric vector of the incident beam 


a, = (sin 20¢) (8 — Og) 
8; = 2a, — Xo (1 — valyo) 
8B. Real part of B,, = 2a,— X»,(1 — v;/79) 
y= By!2C [% xx] 
Yo 
> } Reflection and transmission coefficients 


Transmission coefficient in the absence of reflection. 
= exp (— pD). 


os 
° 


4, DYNAMICAL THEORY OF X-RAY REFLECTION 


(t) 


(K) 


( (Yo) 
i( Yu) 


{ (Yo) 
(vy) 


) (py) 
((by") 


(Do) 
UD) 
ae 
(Dy) 
(n) 
{(k9®) 
(ky?) 
(() 
’) 


ie 


(a/2) 
(2 z/b) 


(y) 


The dynamical theory, in its most general form, seeks to determine the 
interaction of a crystal with an incident beam of radiation and to work out 
the phenomena which accompany this interaction. Here, we are interested 
in the case when the wavelength is of the same order as the lattice spacing 
of the crystal, and in particular when there is only one diffracted ray. The 
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fundamental equations in this case are: 
Ro» — k* 
n° 
R,?— k? 
Re Ps 


a 


Dy a XoDy Xi CD, | 


(3) 


=e, + xem | 


Putting 

Ro = (1+ 59) k and R, = (1 + 5,) k, 
where 5), 5, <1, one obtains the following equation as a condition for the 
compatibility of the two equations in (2): 

45,5,— 2X, (59 + 5,) + C*X,X,— + Xo? = 0. 


We are specially interested in settings close to the Bragg angle. When 
a setting is specified, then there is a definite relation between the incident 
and the diffracted rays, and therefore 5, and 5, are also related. Thus, the 
equation mentioned above becomes a quadratic in 59 alone (or in 4, alone), 
giving two roots for 59. Corresponding to each value of 5), one obtains 
a value of D,/D, for a particular setting and thus there are two wave fields 
inside the crystal, into which the incident beam is split. To calculate the 
transmitted or the diffracted wave outside the crystal, these two wave fields 
have again to be summed up on emergence from the crystal. 


In the application of the dynamical theory to the X-ray reflection from 
a parallel crystal plate, two cases arise according as the diffracted beam 
emerges from the same face as, or the opposite face to, that on which X-rays 
are incident. These two cases are commonly referred to as the “‘ Bragg ” 
case and the ‘‘ Laue” case respectively. However, a different nomenclature 
is called for, since the term “‘ Bragg reflection” occurs often in connection 
with this subject and should not be confused with the reflection in the Bragg 
case. Further, it is also common to use “ Laue reflection ” to mean crystal 
reflection of white X-rays and “ Bragg reflection” for reflection of mono- 
chromatic X-rays. In view of these, the terms “surface reflection’ and 
“internal reflection”? are suggested to refer to the two cases when the 
reflected X-rays come out of the same surface on which the incident X-rays 
strike and when the reflected X-rays pass through the thickness of the crystal, 
respectively. 

If we denote by Dy; and Dy, the magnitude of the electric vectors in the 
transmitted beam corresponding to the two wave-fields 1 and 2, and by 
D,, and Dj, the corresponding quantities in the diffracted beam (of order h) 
then the boundary conditions in the two types of reflection will be for an 
internal reflection 

Dot Dee= Do; Dat Diz= 0 (4) 
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and for a surface reflection. 
Dout+ Du= Do: 

D,, exp [—27/(R,,°Z) D] + D,oexp [—27/ (R,.'z)D] = 0, (5) 
where z is a unit vector along the direction of the thickness. The quantities 
to be determined in each case, are 
for an internal reflection 

Dy exp [— 277 (Ro Z) DJ + Doz exp [— 277 (Ry. 2) D] = Dy“ ) (6) 
D,. exp [— 277 (R,,° Z) D] + Dj2 exp [—- 277 (Ry, z) DJ = D, § 
and for a surface reflection 
D,, exp [— 2ai (Ro,° Z) DJ -- Dos exp [— 27/7 (Ry: 2) D] = Do” L 
Diy a Dy,» = D, j 
5. EXPRESSIONS FOR R AND T IN THE CASE OF INTERNAL REFLECTION 
Laue (1949) has considered this case. The formule (8) and (9) given 
below can readily be derived from his equations. Here, they have been put 
in a form convenient for numerical computation. Further, the rapidly 
varying cosine function in Laue’s formule have been neglected, as these 
only represent interference bands, which arise from interference between 
the two surfaces of the crystal plate. These bands would be very closely 
spaced, the number between two successive orders of Bragg reflection being 
of the order of D/A, which is about 100,000 even with a crystal 0-01 mm. 
thick and even more with thicker crystals. Even a variation of the order 
of A in the thickness of the crystal would obliterate them. It is therefore 
justifiable to neglect these fine variations of intensity and to take an average. 
Eqn. (8) and (9) represent in fact such an average value for the intensity. 


, eXp £() - ‘ ) aD ! 


R = iM Yo YA?) cosh 2D (8) 
Xa 2 cosh? v,, 
l l 
exp | +(~ ¥ _) nD! 
T sn A cosh 2 (nD +- v,) (9) 
2 cosh? vy, 
‘where 
) A 
sinh vy, = 8,/2C [ Vi xx] (10) 
Yo 
and 
7, : ; 
4 ( -) Bye , XarXni7th 
0 : h 0 (1 1) 
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In these equations, 8, is a measure of the departure of the glancing angle @ 
from the Bragg angle 9,, i.e., 8, = 0 corresponds to 4 = 4,. 


It is interesting to note the following deductions which readily follow 
from eqns. (8) and (9): 


4 ‘ 
(a) When |8,| > 2C [% xan | , ie.. for directions sufficiently remote 
0 


rom the Bragg angle, 

R = 0 and T = exp (— pD/y_) = Ty (12) 
as should be the case. This obvious result has been particularly mentioned 
here, because the curves in Laue’s paper may lead to an erroneous impression 
that T increases with 8, on one side of the maximum. His curves | and II 
in Figs. 2 (a), (6), (c) have been drawn without taking into account the 
cosh? v, in the denominator of Eqns. (8) and (9) above. When this is also 
included, the above result (12) follows 


(6) when 8, = 0 


X;, ] l C x; 
T =e R= deer D| (- ed — pe Mt). (13) 
XI diet § Yo a V'y0Yn Xoi 
Thus, for a crystal obeying Friedel’s Law, R = T at the peak of the reflection, 
as has already been pointed out by Laue. 


It is interesting to note that the ratio R/T, depends on the ratio 
X4i/Xoi. i.e., on the ratio of the imaginary components of the structure factors 
of order fh and zero. This ratio would correspond to the F-value, i.e., the 
ratio of the real components, for wavelengths sufficiently far away from 
absorption frequencies. If this is not so, then the peak intensity of an 
internal reflection does not represent the structure factor of the crystal. 


Since X,; < Xp; and C < 1, R cannot exceed the value 4 for a symmetrical 
internal reflection. This is to be expected, because for a non-absorbing 
crystal, R = 4 at the peak. 


~ 


(c) For intermediate values of 8,, it is convenient to introduce the 
quantity y = sinh v,, so that 


1 1 ) y C x, 1 ] 
= i - = zs Ts; «TT | = Ahi I — 14 
’ w [a Yo VMIVIFY  ~ Vra% XG VIF (14) 
Also, we have 
cosh? vy, = (1 + y?) (15) 
From the above definition, ya 8 and hence to (6 — 63), and units of 
y may thus be employed as abscisse in plotting the variation of R and T with 


glancing angle. as has been done by Zachariasen (1945). y= 1 would 
A4 
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correspond to a value of (8 — 4,) equal to the half-width of the correspond. 
ing surface reflection curve if the crystal is non-absorbing (i.c., for the same 
reflection, /). As is well known this is of the order of a few seconds of arc, 
and its magnitude can be readily calculated from @ and the structure factor, 
so that the y-unit forms a convenient unit for representing reflection curves, 


6. CALCULATION OF REFLECTICN AND TRANSMISSION CURVES FOR 
INTERNAL REFLECTIONS 


Laue (1949) has plotted the variation of R and T with angle @ for the 
111 and 200 reflections of rock-salt, with thicknesses of the order of 0-1 mm. 
and using Cu K, radiation. He obtained curves, which exhibit variations 
similar to those observed by Borrmann (1941), namely, an increase and 
decrease of intensity compared to the background on either side of the 
Bragg reflection. The exponent (u D/y) is of the order of 1-0 to 2-0 and is 
of the same order of magnitude as that employed by Borrmann, who used 
quartz crystals 0-2 mm. thick and Fe K, rays. In Campbell’s experiments 
(1951 5) on the other hand, calcite crystals, 0-4 to 1-00 mm. thick were used 
with Cu K,, for which » D/y ~8 to 20. The conditions in the experiments 
of Campbell are not therefore comparable with those assumed in Laue’s 
calculations, and it is not surprising that the latter do not correspond to the 
observations of Campbell. Consequently, the authors have plotted a series of 
reflection and transmission curves, with varying (u D/y), for (101-1) reflection 


of calcite, when the crystal is cut with its surfaces perpendicular to the lattice 
planes concerned. 


For convenience, we write P = » D/y, and a series of curves have been 
plotted for values of P = 0, 0-5, 1-0, 1-5, 2, 3, 4, 6, 8, 20. As mentioned 
in the last section, the abscisse are marked in y-units, y = | actually corres- 
ponding to 3-6” of arc in this case. The relevant data required for the 
calculation have been taken from Zachariasen (1945) and are: 


p = 183 cm? Yo = ¥, = 0-967 
Xo = — 17-44 x 10-8 Xo; = — 0-448 x 10-6 
Xip = —_? 8-57 Xx 10-6 X hi —— 0-429 x 10-6 


The curves in Figs. 1 and 2 correspond only to the normally polarised 
component (C= 1). Since » and y are fixed quantities for the different 
curves, varying values of P correspond to varying thicknesses D, which 
may be readily calculated, since p/y = 189. 


The curve for P =0 is hypothetical, for with an absorbing crystal 
(u #0), P can become zero only if the thickness D vanishes. Actually, 
eqns. (8) and (9) from which the curves have been plotted assumes that the 
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thickness of the crystal is large enough compared with the lattice spacing. 
However, this particular curve may be said to correspond to the case when 

the crystal becomes non-absorbing, other conditions remaining the same. 
All the other curves relate to crystals of varying thickness D. en - 


The following are clear from these curves. Considering first Fig. 1 
(reflection), all the curves are symmetrical about the ordinate for the peak. 
The peak intensity steadily decreases with increasing thickness, but the 
decrease is far less rapid than the corresponding decrease of T,. Thus, for 
P = 0, the peak is 0-5 (T,= 1); this decreases only to about 0-23 for P =2, 
when the transmitted beam has dropped to 0-135. Thereafter the decrease 
in the peak intensity of the reflected beam is much less rapid. In fact for 
P = 20, for which T, will only have a value of the order of 10-%, the peak 
reflection is as large as 0-1, i.e., about a hundred million times the trans- 
mitted intensity for settings away from the Bragg reflection. 














Fic. 1. Reflection curves for varying thickness of the crystal. The parameter P = uD/y. 


The variations in the intensity of the transmitted beam predicted by 
the theory is even more interesting. For P = 0 obviously the transmission 
curve (Fig. 2) is complementary to the reflection curve and exhibits only a 
minimum. However, for small finite values of P (i.e., for small thicknesses), 
it exhibits both a minimum and a maximum and the intensity distribution 
is asymmetric about the peak of the Bragg reflection (i.e., y= 0). On 
approaching the Bragg setting from the left, the transmitted intensity first 
becomes less than the normal value (T,), then it rises and becomes larger 
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Fic. 2. Transmission curves for varying thickness of the crystal. 


than the normal value, finally falling back to Ty. For small values of P, 
i.e., very thin crystals, the decrease predominates, while with increasing 
thickness, the increase becomes larger and larger, until at P = 2 the increase 
is larger than the decrease. For still thicker crystals the minimum is 
imperceptible and only a peak is observed in the transmitted beam. It is 
interesting to note that the peak intensity of the transmitted beam is of 
the same order as that of the reflected beam, although the two peaks do 
not coincide, This is clearer from Fig. 3, which is discussed below, Thus, 
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for P = 20, which would correspond approximately to 1 mm. of calcite, the 
peak value of T is several million times the value of Ty. 


We shall now consider how the theory fits the observations of Campbell. 
For this purpose, Fig. 3 has been drawn to correspond to Fig. 6 of his 











Fic. 3. Reflection (full line) and Transmission (dotted line) Curves of a crystal 

of calcite 0-4 mm. thick. y = 1 corresponds to 3-6’ of arc. 
paper (1951 5), namely, for a crystal of calcite of thickness 0-4mm. It is 
obvious at once that the theory qualitatively explains his observations 
(a) that both the reflected and transmitted beams exhibit sharp maxima 
rising up to several times above the background, (b) that the peak intensities 
of the two are nearly equal to each other and (c) that the peaks are extremely 
sharp. However, quantitative agreement is not to be expected because 
the theoretical curve is for an incident beam of negligibly small divergence, 
while in the experiment, it had a divergence of the order of 15” of arc, which 
is larger than the theoretical width of the reflection or transmission curves. 
Because of this, the observed curves would be much broader than what is 
predicted by theory and the peak intensity would also be lower. When 
these are taken into account, the agreement between theory and experiment 
must be considered good. 


As is clear from Fig. 3, the peaks in the reflected and transmitted beams 
do not coincide. Campbell (1951 5) did observe such an effect, but he 
attributed it to a mis-setting of the crystal. More careful experiments 
should be made to verify whether such an effect is really present or not. 


The curves also show that the peak in the transmitted beam is not 
likely to be observed unless the incident beam has a small divergence, If 
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it has a large divergence, the effect would only be an average over a range of 
angles. This may lead to a net increase or decrease in the transmitted 
intensity over the background according to the thickness of the perfect 
crystal. For very thin crystals, it is a decrease and changes to an increase 
for thicker crystals. This explains the observation of Lonsdale (1947) that 
the deficiency lines observed by her with mosaic crystals became less clear 
and sometimes vanished with perfect crystals. She studied mostly crystals 
of low absorption (diamond and other organic crystals), for which P ~1 
to 2 with a thickness of the order of a few millimetres. Consequently, if 
the whole crystal were perfect, there would be practically no reduction or 
increase in T, while if it is mosaic, a reduction would be observed. Of 
course, the above curves relate to calcite 101-1 reflection and the conditions 
would be different with the crystals studied by Lonsdale (1947), but the 
orders cf magnitude are not likely to be different. The authors have under- 
taken a theoretical study of the variation of R and T with various factors 
such as the structure factor, absorption coefficient, etc., and this would throw 
further light on the subject. 


7. EXPRESSIONS FOR R AND T IN THE CASE OF SURFACE 
REFLECTION 


As mentioned earlier, although Kohler (1933) has considered this case, 
he has made some approximations, e.g., neglecting one of the two wave- 
fields in the crystal, as it has a large attenuation coefficient. It is, however, 
unnecessary to make the approximation and we shall give below a brief 
derivation of the formule, both for the reflection and transmission coeffi- 
cients, following closely the method and nomenclature of Laue (1949). 


In the Bragg case, y,/y, is negative and where necessary this will be 
denoted explicitly by writing it as — | y;/yo|. 
Let 


— 2a, om Xo (1 =o Y1/¥o) 
Analogous to eqns. (29), (30), (31) of Laue (1949), we obtain 








D, —_ xX, Yo $ = 
De (x Vh ye i. i 
where 
ee [4 i 
cosh u = 50 vy, Xb Xi | ‘ (11) 





We have the boundary conditions 
Dat De= Do; 
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and we have to determine the quantities 
D,= Dat Dra; 
Do = Dg: exp [—27i (Roz) D]J4+- Dg, exp [—27/ (Ryg'z) D] (13) 


where z is a unit vector in the direction of the thickness D. If we write 





Ro 2= Ro — kd; 2 Z, (14) 
where 5,, have the values [from Laue, 1949, Eqn. (21)] 
ae | Xo 1. B I [ 2. AC2 YA 7 } 15 
7 “Yo a il eae: Ld sais 
Then, from eqns. (12), we have 
_ exp (u+ x) nin __ Cap — (+ Dine 
0” 2 sinh (u + x) Do"; Da= 2 sinh (u-+ x) Do” (16) 
where 
Substituting in eqns. (13), we obtain 
__ |iX,|| sinhx /? 
ine iX,| |sinh (uw + x)| (18) 
_| . | 5,+ 5.\(? | sinhu |? 
sian — on (C ae ) sinh (u -+ x) (19) 


Unlike the Laue case, it is not possible to put eqns. (18) and (19) in terms 
of real quantities alone. This is so because uw (and x) may be predominantly 
real (i.e., real part S imaginary part) for certain values of y and predomi- 
nantly imaginary for other values.* 

When the reflection is symmetric, i.e., it takes place from lattice planes 
parallel to the surface, and if. Friedel’s law holds (i.e., |x,| = |x,}), then 
eqns. (18) and (19) reduce to 


) 


sinh x (|? | sinhu (|? 
R= sinh (u + x) ia ‘sinh (u + x) (20) 
These niay be compared with the corresponding expressions obtained by 
one of the authors (Ramachandran, 1942, 1944) for a finite parallel-sided 
crystal plate having 7 lattice planes parallel to its surface. The latter are 

_ | sinh nf 4 | 

R = |aah (e + mp) T= (neta (21) 
It is clear that there is an exact correspondence between the two, and an 
examination of the significance of the quanitties a and nf in eqns. (21) 


——————— 
* e.g., when there is no absorption, u is pure imaginary for —1< y< +1 and 
completely real for values of y outside these limits. 


sinh a (2 
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shows that they are the same as u and x. This again confirms the statement 
made by Ramachandran (1948) that the Ewald and Darwin methods lead to 
identical results, in spite of their apparently very different mathematical 
approaches. 

SUMMARY 


The paper deals with the theory of an interesting phenomenon (which 
has been designated as “anti-reflection ’’) that the intensity of the trans. 
mitted beam may exhibit a peak larger than the background when a Bragg 
reflection occurs in an absorbing crystal. The theory is based on the 
dynamical theory of Ewald and Laue. It comes out that the effect is due 
to a decrease in the effective absorption coefficient of the crystal near the 
Bragg reflection, and to the consequent increase in the transmitted intensity 
predominating over the loss of energy by reflection. The anti-reflection 
peak becomes more prominent, the greater the thickness of the crystal. The 
results of the theory are found to be in accord with the previous observations 
of Borrmann and Campbell. The theory further predicts that the peaks 
in the reflected and transmitted beams would not be coincident and this 
requires further verification. 
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